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1.0 INTRODUCTION 

This report summarizes the results obtained on the NASA-ARMY Add-on 
to the "Preparation of High Purity Low Dislocation Density GaAs Single 
Crystals" program (NAS3-22224), The primary goal of the add-on program was to 
increase the scope of the original program to determine the materials syn- 
thesis and growth conditions necessary to produce high purity low dislocation 
GaAs by the liquid encapsulated Czochralski (LEC) technique. Attention was 
focussed on two principal areas. First, the initial study of crystal growth 
parameters that control the dislocation density in large-diameter (3-inch) 
crystals was expanded. Low dislocation densities are particularly important 
for solar cell applications requiring high performance and radiation hard- 
ness. Second, the predominant chemical impurities and electrically active 
centers in undoped LEC material and their influence on the electrical compen- 
sation were evaluated. 

Three important accomplishments have been achieved in our laboratory 
as a result of research conducted through the NASA-Army and our IR&D pro- 
grams. First, and most important, a crystal growth process was developed for 
producing semi-insulating CaAs with a demonstrated high yield of device- 
quality material. This process involves the control of the melt stoichiometry 
on the As side of the stoichiometric composition. Evaluation of the material 
grown by this process showed that the semi-insulating properties are uniform, 
and thermally stable during device processing. LEC GaAs wafers have been 
successfully used at Rockwel 1 -MRDC to fabricate integrated circuit devices 
such as RAMS, multipliers, and synthesizers. The implementation of this 
crystal growth process and the use of the resulting GaAs has had considerable 
impact on the development of GaAs integrated circuit technology. Rockwell - 
MRDC is the first laboratory to upgrade the device fabrication line to process 
3-inch GaAs wafers as a result of this work. 

Second, studies of the structural perfection of LEC GaAs crystals in 
this program have led to an additional 10% reduction in the dislocation 
density compared to the results achieved during the first phase of the pro- 
gram. Reduced dislocation densities were found to be favored by the use of 
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thick B 2 O 3 layers* low ambient pressures, and controlled melt stoichiometry. 

By controlling the appropriate crystal growth parameters, the most struc- 
turally perfect 3-inch-diameter LEG GaAs crystals yet reported were produced 
in this laboratory. 

Finally, it was found that the control over the melt stoichiometry is 
also essential in reducing the frequency of twinning and polycrystallinity 
during crystal growth. This understanding dramatically increased the yield of 
single (100) GaAs crystals required for device fabrication. 

Six crystals were grown as a part of the NASA-Army program (Fig. 1). 

Two of these crystals (R30 and R39) were delivered as whole single crystals to 

Mr. Thomas AuCoin in fulfillment of the contract obligations. The four 
remaining crystals (R21, R24, R33, and R34) were sliced into wafers. Four 
polished wafers each from the front and tail of each of these four crystals 
were also delivered to Mr. AuCoin. Other wafers were characterized using 
various destructive and nondestructive techniques such as secondary ion mass 
spectrometry, Hall effect measurements, and preferential etching. In 
addition, a number of these semi-insulating substrates were processed into IC 
devices for the experimental evaluation of cause-effect relationships between 
substrate properties, device processing, and device performance. 

The organization of the report is as follows. In the second section, 
the materials synthesis and growth techniques used to produce the crystals are 
described, with emphasis on new procedures used and new knowledge gained since 
the original program. Various techniques used to characterize the material 
are described in the third section. The central results concerning the 
compensation mechanism and dislocation studies, are described in the fourth, 
and fifth sections, respectively. The sixth section describes the use of 
these semi-insulating materials in device fabrication. 
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2.0 MATERIALS SYNTHESIS AND GROWTH TECHNIQUES 

In this section, the processes Involved In materials synthesis and 
growth of LEG GaAs are desc>*1bed. 

2.1 Growth Configuration 

All crystals were grown in Rockwell International's Melbourn high- 
pressure LEC system at the Thousand Oaks MRDC laboratory. The configuration 
of the LEC system, shown schematically in Fig. 2, consisted of a GaAs melt 
contained In a high purity pyrolytic boron nitride (Union Carbide) crucible. 
The boric oxide (B2O3) encapsulant floats on the top surface of the melt. In 
addition, a thin film of B2O3 coats the entire surface of the crucible due to 
the high temperature wetting characteristics of this material. The ^2^3 also 
wets the growing crystal. Thus, the GaAs melt is completely sealed, suppress- 
ing As evaporation and shielding .the melt against contamination from the cru- 
cible and the growth ambient. 

2.2 Growth Process 

The major steps In the crystal growth operation included loading of 
the charge, heat up, synthesis, equilibration, seeding, necking, cone growth, 
and pulling of the full diameter ingot. The crucible was loaded with approx- 
imately 1400 g 6-9's Ga (Ingal International), 1500 g 6-9's As (Cominco), and 
a 500 g pre-formed B2O3 disk (Puratronic) with a known moisture content, as 
shown in Fig. 3. Solid Ga, which is solid to just above room temperature, was 
loaded on top of the As so that the liquid Ga served to encapsulate the As. 
Starting with a chamber pressure of 600 psi, the crucible was heated to 
between 450 and 500®C, at which point the B2O3 melted, flowed over the charge 
of Ga and As, and sealed at the crucible wall. The synthesis reaction 
(Ga^^qy^jj + ASjjq^^^I = GaAs^^^^^) occurred at about 800®C. The presence of the 
B2O3 and the use of high overpressures (^1000 psi) prevented significant loss 
of As due to sublimi nation and evaporation during and subsequent to synthe- 
sis. The melt was then equilibrated at the starting temperature and the 
growth procedure begun. 
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Fig. 2 Cross section of the crucible for the LEC Growth system, 
showing the location of the B^Oj during growth. 
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Growth was initiated by dipping the seed, which was held on the pull 
shaft, through the B2O3 and into the melt. The crystal was grown by gradually 
withdrawing the seed from the melt* The diameter was gradually and control- 
lably increased to full dimension. The seed and the crucible were rotated in 
the same direction at 6 and 15 rpm, respectively. 

The growth parameters and physical properties of the NASA-Army crys- 
tals are given in Table 1. 


Table 1. Growth Parameters and Physical 
Properties of NASA-Army Crystals 


Crystal 

No. 

Crucible 

Material 

B2O3 Water 
Content, ppm 

Growth 

Direction 

Weight 

(kg) 

Diameter 

(mm) 

R21 

PBN 

495 

<100> 

2.2 

71.1 i 7.1 

R24 

M 

170 

II 

2.3 

82,4 t 5.5 

R30 

II 

170 

II 

2.4 

81.3 t 3.0 

R33 

M 

160 

II 

2.1 

See text 

R34 

11 

' 304 

M 

2.4 

87.5 + 6.3 

R39 

II 

335 

II 

2.5 

84.5 +2.6 


As will be discussed later in this report, control over the melt 
stoichiometry is essential to achieve reproducible growth of undoped semi- 
insulating LEC GaAs. An important element of the pre-growth process is method 
by which the charge of elemental Ga and As is heated from room temperature to 
the melting point of GaAs, In the earlier experiments with PBN crucibles it 
was found that significantly more As was lost from PBN crucibles compared to 
qua rt? crucibles during the heat up cycle. The typical weight loss from a 
3000 gram charge contained in a quartz crucible ranged from 10 to 20 grams, or 
less than 1%. The composition of the charge could effectively be controlled 
under these circumstances. On the other hand, the As weight loss from PBN 
crucibles could be Several times the loss factor of quartz crucibles, leading 
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to significant variations in melt stoichiometry as a crystal was pulled from 
the melt. (The weight loss was determined by comparing the weight of the 
initial charge to the weight of the crystal and the charge remaining in the 
crucible after growth.) 

Initially, the difference in behavior between the two crucibles was 
believed to be a result of the different heat-transfer characteristics of the 
two materials; the higher heat transfer of the PBN at Tower temperatures 
(<600*C) would lead to premature heating of the As before the 82 O 3 had com- 
pletely melted. For this reason, a number of ixperiments were conducted at 
the beginning of this program in which the heating rate was reduced to allow 
sufficient time for the B 2 O 3 to melt. The dependence of the weight loss on 
heating rate in experiments 21, 24 and 30 are given in Table 2. The results 
show that the weight loss Increases as the heating rate decreases, which is 
opposite to the expected behavior based on the original model for the weight 
loss. Therefore, it was concluded that the difference between the crucibles 
is vyf; the heat-transfer characteristics, but rather, the wetting character- 
l||?C:s of the B 2 O 3 . Since quartz is an (ionic) oxide, whereas PBN is a 
(covalent) semiconductor, different wetting characteristics for B 2 O 3 the two 
materials are to be expected. The weight loss decreases as the heating rate 
increases because the wetting to the PBN is Improved at higher temperatures, 
allowing the B 2 O 3 to seal before significant amounts of As are lost. Using 
high heating rates in combination with additions of either excess As or Ga in 
the initial charge, we can now routinely produce near-stoichiometric 3 kg 
melts with the desired composition offset (either As-rich or Ga-rich) with 
respect to the stoichiometric composition. 
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Table 2. 

Dependence of As Weight Loss 
Rate During Heat-up cycle 

on Heating 

Crystal No. 

Average Heating 
Rate, (®C hr"*) 

Weight Loss 
(grams) 

R21 

14S 

140 

R24 

109 

232 

R30 

800 

91 
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3.0 CHARACTERIZATION TECHNIQUES 


The techniques used to characterize the crystalline perfection* chetn- 
ical impurities* and electrically active centers In the LEC GaAs are sum- 
marized in Table 3. 


Table 3. Summary of Techniques Used to Characterize LEC GaAs 


Material 

Parameter 

Characterization 

Technique 

Consents 

Dislocation 

Density 

Preferential Etching 

1) For (100) wafers 

2) 25 min in KOH at 400’C 

Chemical 

Impurities 

Secondary Ion Mass 
Spectrometry (SIMS) 

1) Charles Evans and Assoc. 
San Mateo, CA 

2) Good for metals and donor 
impurities in GaAs 

Carbon 

Localized Vibrational 
Mode Far- Infra red 
Absorption (LVM) 

1) Measurements made at 77K 

2) Reference - M. R. Brozel , 
J. B. Clegg, and 

R. C. Newman (1978) 

J. Phys. DJ^, 1331 

Deep Levels 

Photoluminescence 

Made by Phil Yu, Wright State 
University, Dayton, Ohio 

EL2 (deep donor) 

Optical Absorption in 
1.0-1.4 urn range 

Reference - G. M. Martin 
(1981), Appl. Phys. Lett. 39, 
747 

Resistivity, 
Free Carrier 
Concentration, 
Hall Mobility 

Hall Effect Measurements 

1) 300K 

2) Dark 

3) AuGerNi Contacts 

Activation 
Energy of 
Electron 
Concentration 

High-temperature 
Hall measurements 

300-420K 
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4.0 COMPENSATION MECHANISM OF SEMI -INSULATING LEC GaAs 

Results from the NASA-Army and Rockwell IRffl programs indicated that 
the key to the reproducible growth of undoped semi -insulating GaAs by the 
liquid encapsulated Czochralski technique is the control over the melt 
stoichiometry. Studies in the crystal growth IR&D program revealed that the 
electrical compensation in the semi -insulating material is controlled by the 
balance between EL2 deep donors and carbon acceptors. Furthermore* the 
incorporation of EL2 depends on the melt stoichiometry, increasing as the As 
atom fraction in the melt increases. As a result, semi -insulating material is 
obtained only above a critical As fraction in the melt. Below the critical 
composition, the material is p-type. The results of this work, including the 
identification of the predominant chemical impurities and defects, the 
development of the compensation model, and the evaluation of the material in 
terms of the model, are discussed in detail in the two publications given in 
Appendices I and II. A summary of these results is given in the following 
section. 


4. 1 Importance of Melt Stoichiometry 

The resistivity of the LEC GaAs is controlled by the melt stoichiom- 
etry {Fig. 4). The material is semi -insulating above, and p-type below, a 
critical melt composition. For example, referring to the electrical proper- 
ties of the NASA-Army crystals in Table 4, and the melt stoichiometry corres- 
ponding to the front and tail of these crystals in Table 5, it can be seen 
that the resistivity of the crystals grown from As-rich melts (R33 and R34) 
are semi-insulating from front to tail. However, crystals grown from Ga-rich 
melts either undergo a transition from semi -insulating to p-type at some point 
along the crystal, as was the case with R21 and R24, or are p-type throughout 
(Fig. 5). These results show that the transition to p-type conduction is 
controlled by the melt stoichiometry and is not related to the normal segre- 
gation of a background acceptor impurity toward the tail of the crystal. 
Otherwise, the tail of the As-rich-grown crystals would be p-type as well. 
Detailed resistivity profiles of R21 and R24, shown in Table 6, pinpoint the 
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042 0,44 0 4 6 048 0.50 0.52 0.54 

ARSENICATOM FRACTION IN MELT 


Fig, 4 Dependence of electrical resistivity of LEG GaAs on melt 

stoichiometry. Semi -insulating material is obtained abovei, 
and p-type (low resistivity) material below a critical melt 
composition of about 0,475 atom fraction. 
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Table 4. Electrical Characteristics of NASA-ARMY Crystals 


Sample Ho.* 

Resistivity 

a-cm 

Type 

Free Carrier 
Concentration 
cm"* 

Hall 

Sheet 

Resistivity 

n/O 

R21T 

1.5 X loB 

n 

2.1 X 107 

1962 

2.5 X IQ9 

R21F 

0.60 

P 

3.9 X 1016 

278 

- 

R24F 

1.3 X 107 

n 

O 

f-4 

X 

00 

2774 

2.5 X 109 

R24T 

0.65 

P 

2.9 X 1016 

328 

- 

R30 

Delivered whole 




R33F 

4.2 X 1q7 

n 

3.4 X 107 

4433 

1.3 X 109 

R33T 

7.3 X 106 

n 

1.9 X 108 

4614 

1.0 X 10® 

R34f' 

6.1 X 107 

n 

2.5 X 107 

4189 

1.0 X 109 

R34T 

3.4 X 107 

n 

3.8 X 107 

4812 

• 

X 

o 

00 

R39 

Delivered whole 





* F * Front of Ingot 
T * Tall of Ingot 
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Fig. 5 Resistivity profiles for LEC GaAs crystals grown from Ga-rich 
and As-rich melts. 
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critical melt compositidn of the semi-insulating-to-p-type transition at 0.475 
As atom fraction in our material. 

4«2 Chemical Impurities and Electrically Active Centers 

To gain a better understanding of the compensation mechanism in the 
semi-insulating material, the predominant chemical Impurities and electrically 
active centers were identified. The concentration of background impurities in 
the NASA-Army crystals, shown in Table 7, are consistent with SIMS and LVM 
studies of approximately 20 crystals grown in this crystal puller. Boron in 
the predominant chemical impurity in the material (6 x 10^® to 2 x cm"^). 
However, boron is isoelectronic with Ga and is not electrically active in this 
concentration range. Carbon is the predominant electrically active chemical 
impurity (3 x 10^® to 1.5 x cm"^). Note that the Si concentration is 
typically below the I x lo* cm" level, more than a factor of about 10 lower 
than in Bridgman material analyzed. 

The predominant deep center in the semi-insulating material, studied 
as a part of the IR&D work, is the deep donor commonly referred to as EL2. 

This center was identified by optical absorption in the 1.0 to 1.4 urn range 
(Fig. 6) and by the activation energy of the electron concentration (Fig. 7), 
which was found to be 0.75 t 0.02 V. For the first time, it was shown that the 
EL2 concentration in the bulk material depends on the composition of the melt 
(Fig. 8). 
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Fig. 6 Typical optical absorption spectrum of GaAs containing EL2. 
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Fitj. 7 Temperature-corrected electron concentration - temperature 

curve of undoped semi-insulating 6aAs giving activation enfergy 
of the free carriers. 


20 



CONCENTRATION OF EL2. cm*^ 








Rockwell International 


25x10 


15 


15 


" 20 X 10 


15x10^® 


10x10 


15 


5x 10 


15 


T r 


T r 


STOICHIOMETRIC 

COMPOSITION 


$ 


MRDCI1^14021 
“—1 


0.47 0.48 0.49 0.50 0.51 0.52 0.53 0.54 

ARSENIC ATOM FRACTION IN MELT 



I 

I 


Fig. 8 Dependence of EL2 concent^'ation as determined by optical 

absorption on the melt stc^lchiometry. The concentration of 
EL2 increases from 5 x 10^5 cm“^ to 1.7 x 10^^ cm*^ as the 
As atom fraction increases from about 0.48 to 0.51, and 
appears to saturate as the As concentration increases further 
to 0.53* 
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5.0 REDUCTION OF THE DISLOCATION OENSITT 

The effects of several crystal growth parameters on the dislocation 
density were evaluated In the original NASA program. It was shown that 
reduced dislocation densities are favored by employing low dislocation seeds 
and Dash-type seed necking, by maintaining good diameter control, and by 
growing the crystal with a cone angle greater than about 25“. During the 
course of the NASA-Army Add-on program, the additional effects of ambient 
pressure, melt stoichiometry, and 8303 height were studied. More detailed 
investigations of the effect of the cone angle were also made. The results 
are discussed In the paper given In Appendix III, The pertinent results are 
summarized In this section, and the experimentally determined dislocation 
densities are shown In Table 8 . 

5,1 Cone Angle 

In the previous program. It was reported that the dislocation density 
Is virtually Independent of the cone angle 0 for 20® < 9 < 70®, and that 
unusually high densities occur In crystals grown with a "flat top" ( 0 ® < 9 < 20 ®), 
During the Add-on program, R34 was grown with a 0 ® cone angle (which we had 
never done before), representing the Ideal "flat-topped" crystal. This 
experiment was critical to our unders^ nding of the origins of dislocations in 
LEC GaAs. 

The etch pit density (EPD) In the front of R34 was unusually high 
(low 10® cm"^ range) compared to the other crystals. In addition, the longi- 
tudinal EPD (along the crystal) In the front of the crystal was inverted, 
first decreasing, and then increasing toward the tail. Also, during crystal 
growth, the crystal began to expand rapidly when the top of the crystal 
emerged from the B 2 O 3 encapsulating layer, leaving a bulge at the distance 
from the front of the crystal equal to the height of the 8303 . 

These observations Indicate tnat the top of the crystal was subjected 
to significant additional cooling as the crystal emerged from the B 3 O 3 
layer. As a result, the flat-top crystal experienced "thermal shock" as it 
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Table 8, Dislocation Density and Distribution 
In NASA-Army Crystals 





Dislocation Density, 

cni 



Sample No. 


Front 




Tall 


1 

2 

3 

Region 

1 

2 

3 

R21 

2.3e4 

3.3e4 

3.0e5 


l.OeS 

l.OeS 

2.4e5 

R24 

l.2e4 

2.6e4 

l.leS 


♦ 



R30 




Deliverable 




R33 

6e3 

1.8e4 

9.6e4 




** 

R34*** 

l.leS 

2.4e5 

2.7e5 


1.3e5 

2.3e5 

1.6e5 

R39 




Deliverable 





* Data not obtainable because the crystal twinned. The orientation of the 
wafers obtained from the tall of the crystal was not (100). Our preferen 
tial etch can be used only for (100) orientations. 

** Qa Inclusions, which formed as a result of surface degradation under low 
pressure, led to extremely high EPD in tall uncharacteristic of as-grown 
material (see text). 

*** Unusally high EPD In front resulted from flat-top growth (see text)* 
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suddenly emerged from the B2O3. Presumably, the inverted dislocation density 
along the crystal resulted from the formation of a “secondary” distribution of 
dislocations, which formed after, and added to the grown-in distribution. 

Detailed studies of longitudinal cross sections of crystal cones 
varying in cone angle from 20 to 65® further support the conclusion that addi- 
tional dislocations form near the top of the crystal as the crystal emerges 
from the B2O3. 

The radial EPO distribution (across wafers) is “W'-shaped (Fig. 9). 

In addition, the EPD at a particular position in the crystal body at a 
constant distance from the center increases toward the tail. However, the EPD 
In the cone region actually increases, and then decreases before the crystal 
reaches full diameter, resulting in a local maximum of the EPD. The EPO value 
at the maximum increases, and the position of the maximum is closer to the 
seed, as the cone angle decreases. This behavior is consistent with the 
results obtained on the flat-top crystal, R34. Indeed, R34 represents the 
limiting case of a 0® cone, in which the EPD value at the maximum is the 
highest observed in any of these crystals, and the position of the maximum 
within the crystal lies directly below the seed. 

These results are important because they show that the convective 
heat transfer from the surface of the crystal to the high pressure ^r ambient 
is greater than the heat transfer to the liquid B2O3 encapsulating layer. As 
a result, radial temperature-gradient-induced stress increases as the crystal 
emerges from the B2O3. Since the crystal above the encapsulating layer is hot 
enough to undergo plastic deformation, additional dislocations are created. 

Two additional growth parameters were varied to reduce the heat 
Tosses through the ambient gas, the B2O3 layer, and the ambient pressure. 

5.2 B2O3 Thickness 

The EPD in the front of the crystals was found to decrease as the 
thickness of the B2O3 increases as shown in Table 9. In view of the discus- 
sion above concerning the relative heat loss to the ambient compared to the 
heat loss to the B2O3, the reduced EPDs brought about using thicker B2O3 
layers probably 


24 

G/3748A/cb 


ETCH PITS DEMSITY |cm^ 






Rockwell International 


ORIGINAL PAGE IS 
OF POOR QUALITY 

Table 9. Effect of B 2 O 3 Height on Dislocation Density 


Ingot No. 

Weight of*" 
B 2 O 3 


Front 

EPO (cm“2) 

Tall 




1 

to 

• 

X 

10^ 

8.0 X 

10^* 

R32 

270 gm 

2 

1.0 X 

10® 

1.1 X 

10®* 



3 

1.5 X 

10® 

2.0 X 

10®* 



1 

1.3 X 

10^ 

1.4 X 

10® 

R35** 

390 gm 

2 

2.8 X 

10^ 

2.2 X 

10® 


• 

3 

1.7 X 

10® 

2.4 X 

10 ® 



1 

1.2 X 

10^ 

9.0 X 

104 

R29** 

500 gm 

2 

1.7 X 

10^ 

1.0 X 

10 ® 



3 

2.5 X 

10® 

2.2 X 

10® 


+: Other growth parameters are similar. 

*: '3/4 Ingot length area. 

Ingots grown under Rockwell IR4D program. 
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result from more effective thermal Isolation between the hot region of the 
crystal and the Ar ambient due to an increase in the effective ambient 
temperature. These results suggest that further increases in B 2 O 3 layer 
thickness will lead to even lower EPDs. 

5,3 Ambient Pressure 

Crystal R33 was grown at an ambient pressure of 50 psi, which is low 
compared to the standard pressure of 300 psi. The EPO at the front of this 
crystal was compared to that of another crystal grown under very similar con- 
ditions, except that the ambient pressure was 300 psi. The results, shown in 
Table IV of Appendix III, show an average reduction in EPO of 45% in the 
crystal grown at low pressure. This behavior can be explained. 

The heat transfer coefficient of the crystal -ambient surface is 
expected to increase as the square root of the ambient pressure. Therefore, 
we could expect a reduction in the heat transfer coefficient by a factor of 
about 2.5 by reducing the pressure from 300 to 50 psi. The experimental find- 
ing of a 50% reduction in EPO is consistent with the theoretical prediction. 

Although further reductions in ambient pressure would lead to 
corresponding reductions in EPO, serious problems could arise due to thermal 
degradation (preferential As evaporation) of the surface of the crystal. In 
fact, significant surface degradation was observed on R33, Since As preferen- 
tially evaporates, liquid Ga remains on the surface which can migrate into the 
crystal yielding poor-quality crystals. Therefore, the ambient pressure 
cannot be reduced to too low a level in an effort to achieve low dislocation 
material without sacrificing the entire crystal. Further work is needed to 
identify the minimum acceptable ambient pressure in this system. 
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6.0 REDUCED INCIDENCE OF TWINNING 

A major problem which had affected the yield of GaAs material suit- 
able for device processing was the incidence of the twin formation. Twinning 
causes changes in the crystallographic orientation of the material and can 
lead to polycrystallinity and grain-boundary formation. 

We have found that control over the melt stoichiometry is important 
to prevent twin formation, as described in Appendix IV. Twenty crystals were 
grown from stoichiometric and nonstoichiometric melts. The incidence of twin- 
ning was shown to be reduced significantly when crystals are grown under As- 
rich, near-stoichiometry conditions as Indicated in Table 10. We are 
presently experiencing better than a 90% yield of single-crystalline material 
using this recipe. 
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7.0 APPLICATION OF LEC GaAs IN FABRICATION 

Several wafers obtained from the NASA-Army crystals were used in 
digital integrated circuit processing with excellent results. Both discrete 
devices such as FETs, as well as circuits, such as 9-stage ring oscillators 
and RAMs (random access memories) were fabricated, A micrograph of a RAM 
fabricated under our DARPA-supported RAM development program is shown in 
Fig. 10. We have also evaluated the uniformity of implanted layers activated 
directly in the material. Good uniformity is important to fabricate high 
performance circuits. 

The uniformity of implanted layers in the NASA-Army material was 
evaluated by mapping the depletion-voltage of Schottky barriers across 
quarters of 3-inch wafers by C-V profiling, as shown in Fig. lla-d. (The 
depletion voltage is analogous to the pinchoff voltage in an FET.) The 
average depletion voltage V, standard deviation of the depletion voltage Oy, 
and the number of diodes tested are summarized in Table 11. We have also 
fabricated arrays of 1-um-gate FETs in our process line on several wafers 
obtained from R34. The standard chip size for processing is 1-inch squares. 
Maps pf the actual pinchoff voltage Vp of four of these wafers are shown in 
Fig. 12a-d, and values of V, Oy, and n are given in Table 12. 

Our experience at Rockwell in the GaAs IC technology indicates that 
the uniformity of Vp (and V^j) in processed wafers reflects both the material 
and the processing. We typically see values of Oy of about 50-70 meV in 
processed circuits, which is an acceptable tolerance for depletion-mode 
devices. An exceptionally low value of 25 meV was obtained on one of the R34 
wafers CR34A (Fig. 11a, Table 12). In fact, to our knowledge, this is the 
best result ever achieved at any laboratory by direct ion implantation. This 
result illustrates the degree of uniformity achievable by direct implantation 
in undoped LEC GaAs. Although all four values of Cy for R34 are excellent, 
the data could indicate a significant contribution of the device processing to 
the variability of Oy in the circuits because these four wafers were processed 
together. The relatively high value of the standard deviation of the deple- 
tion voltage on the quadrant map of R34 (Fig. lid) would also seem to be 

I 
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Fig. 10 Photograph of RAM processed on LEG GaAs. 
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Fig. 11 (a-d) Depletion voltage maps of NASA-Army crystals for 320 keV Se 

Implanted through 750A SljH;. to a dose of 2.5 E12/cm^. 
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Fig. 12 Pinchoff voltage maps of 
crystal R34. 
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Table 11. Summary of Studies of Depletion Voltage 
Uniformity Across 3-lnch Wafers (C-V) 


Crystal 

Noc 

Average 
Depletion 
Voltage, V 
7 

Standard Deviation 
of Depletion 
Voltage, mV 

Ov 

Number of 
Diodes 
Tested 

n 

R21 

1.360 

132 (9.7X) 

145 

R24 

1.644 

120 (7.3%) 

134 

R33 

0.920 

47.7 (5.2%) 

170 

R34 

1.430 

109.9 (7.7%) 

175 


♦The variations of 7 from sample to sample are due primarily to the fact that 
special processing was required to make these measurements. The top surface 
of the Implanted wafers was Inadvertently etched when the SI3N4 cap was 
removed, leading to variations in the thickness of the Implanted layer and 
corresponding variations In 7. This etching step could also have influenced 
the uniformity across each wafer. 


Table 12. Pinchoff Voltage of FET Arrays on Simultaneously 
Processed 1-Inch Wafers from NASA-Army Crystal R34 


Wafer 

No. 

Pinchoff Voltage (Vq) 
V ^ 

Standard Deviation 

of *p (o,) 

mlV 

Number of 
FETS, 

Operating/Total 

224 

0.9441 


24.9 (2.6%) 

71/72 

223 

1.131 


60.0 (5.3%) 

72/72 

222 

1.040 


41.0 (3.9%) 

70/72 

221 

1.015 


47.2 (4.6%) 

72/72 
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inconsistent with the 25 meV standard deviation of Vp of the FET arrays fabri- 
cated from the wafers from the same crystal. The apparent inconsistency could 
also be an indication of process-induced variations. 

The wafer quadrants used for Vj mapping studies (Fig. 11) could not 
be processed simultaneously because of their size and the fact that special 
processing was required. As a result, it is difficult to draw hard conclu- 
sions concerning possible substrate-related variations. Values of are 
higher near the edge of all four samples in the <100> direction compared to 
the center of the wafer. This trend correlates with the general trend of the 
dislocation density, the dislocation density also being higher near the edge 
in the <100> direction. However, the opposite behavior has been observed in 
other quadrant mapping studies. We are now proceeding to relate the uniform- 
ity of device parameters across wafers to the properties of the material. 

Finally, we note that very high yields of operating devices has been 
achieved using R34. The yields of operating FETs experienced in the pinchoff 
voltage mapping studies exceeded 97% (see Table 13). We have also fabricated 
five types of 9-stage ring oscillator circuits on R34 for the purpose of opti- 
mizing the circuit design. The yield exceeded 92% on each wafer, and the net 
yield (170/180 operational) was 94%, which is excellent. Although the yield 
factor may generally be more of an indication of "good" processing, these 
results do suggest that direct material-related device failure is probably not 
an important consideration at the present stage of development of ICs. 


Table 13. Yield of Five 9-Stage Ring Oscillators 
Fabricated on NASA-Army Crystal R34 


Circuit No. 

Devices 

Operating/total 

Yield, % 

1 

35/36 

97 

2 

35/36 

97 

3 

34/36 

94 

4 

33/36 

92 

5 

33/36 

92 
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8.0 SUMMARY 

The add-on Army program resulted in substantial progress in the area 
of LEG growth of high purity, low dislocation undoped GaAs. The results from 
this program, combined with other findings, contributed to the development of 
a compensation model for undoped, semi-insulating GaAs. By investigating the 
effects of cone angle, boric oxide encapsulant thickness, and ambient pres- 
sure, we developed procedures for low dislocation density large diameter LEG 
GaAs. In addition, the incidence of twin formation was reduced by utilizing 
As-rich melts. 

These developments in the LEG GaAs growth process have resulted in 
substrates with superior electrical and crystalline properties. The outstand- 
ing device yield and uniformity data on these substrates reflect these recent 
Improvements. 
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Stolchiometf7-controlled compensation In liquid encapsulated Czochralski 

QaAs 

D. E. Holmes, R. T. Chen, K. R. Elliott, end C. Q. Kirkpatrick 

Rockwell International Microelectronics Research and Development Center, Thousand Oaks, California 
91360 

(Received 13 August 1981; accepted for publication 9 October 1981) 

We show that the electrical compensation of undoped GaAs grown by the liquid encapsulated 
Czochralski technique is controlled by the melt stoichiometry. The concentration of the deep 
don'ir EL2 in the crystal depends on the As concentration in the melt, increasing from about 
Sx iO'^ cm~^ to 1.7 X 10'^' cm~^ as the As atom fraction increases from 0.48 to O.Sl. 

Furthermore, we show that the free-carrier concentration of semMnsuIating GaAs is determined 
by the relative concentrations of EL2 and carbon acceptors. As a result, lemMnsulating material 
can be obtained only above a critical As concentration (0.475*atoffl fraction in our material) where 
the concentration of EL2 is sufficient to compensate residual acceptors. Below the critical As 
concentration the material is p type due to excess acceptors. 

PACS numbers: 72.80.Ey, 81.10.Fq, 78.50.Ge 


Undoped semi*insulating GaAs grown by the liquid en- 
capsulated Czochralski (LEG) technique is currently receiv- 
ing attention for applications in integrated circuits'*^ fabri- 
cated by direct ion implantation. An understanding of the 
electrical compensation mechanism is necessary to allow 
controlled and reproducible crystal growth. In the present 
investigation we have analyzed the electrical and optical 
properties of 1 2 LEG GaAs crystals grown from stoichiome- 
tric and nonstoichiometric melts. The results show that the 
stoichiometry of the melt controls the electrical compensa- 
tion of the crystal through incorporation of EL2, a defect 
that has been implicated'^ as the compensation-controlling 
center. 

The GaAs crystals used in the investigation were grown 
in the <100) direction from pyrolytic boron nitride (PEN) 
crucibles using a high pressure LEG crystal puller, The crys- 
tal diameter and weight ranged from 2-3 in. and 2.0-2.S kg, 
respectively. Dry (<500-ppm HjO) B2O3 encapsulant (John- 
son-Matthey) was used for all growth experiments. The melt 
was prepared by in-situ synthesis, to achieve high purity, by 
starting with a charge of 6-9’s As (Gominco) and 6-9's Ga 
(Kawecki). To make an accurate determination of the initial 
melt composition, it was necessary to take into account the 
loss of As from the charge during the heat-up cycle. We 
found that the As loss resulted from incomplete wetting of 
the BjOj to the PBN crucible before synthesis. This weight 
loss was determined by comparing the weight of the initial 
charge with the weight of the crystal and the charge remain- 
ing in the crucible after the growth process. The As concen- 
tration of the initial melt was effectively varied from 0,46 to 
0.51-atom fraction. Samples for characterization were ob- 
tained along the length of each crystal. The melt composi- 
tion for each of these samples was determined by adjusting 
the initial melt composition for the crystal weight at the time 
of growth. The crystal weight and length were recorded dur- 
ing growth as a function of time. 

The material was characterized by secondary ion mass 
spectrometry'* (SIMS), localized vibrational mode (LVM) 
far-infrared^ absorption, room temperature and variable- 


temperature (77-400 K) Van der Pauw Hall effect measure- 
ments, near-infrared optical absorption, and variable tem- 
perature photoconductivity. The predominant background 
impurity as determined by SIMS was boron (1 X lO'^cm"’’- 
2 X 10” cm"*) and by LVM carbon (3 X 10” cm"-*- 
1.5X 10” cm"*). The concentrations of other background 
impurities were typically Si < 2 X 10” cm"*, S < 2 X 10'* 
cm"*, Fe < 3 X 10'* cm"*, and Mn < 1 X 10'* cm"*. 

The electrical resistivity of the material was found to be 
a strong function of the melt stoichiometry, as shown in Fig. 
1 . These data were obtained from the front and tail of each of 
the 12 crystals and from detailed electrical profiles of 5 crys- 
tals. Semi-insulating material was obtained above, and p~ 
type material below a critical As concentration of 0.475- 
atom fraction. The free-hole concentration and Hall mobil- 
ity of the material grown from Ga-rich (As concentrations 
below 0.50) melts below the critical composition are in the 
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FIG. 1. Dependence of electrical resistivity of LEG GaAs on melt composi- 
tion, Semi-insulating material is obtained above, andp-type material below 
a critical melt composition of 0.<475 As atom fraction. 
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rtnfeofl-3XlO'‘cm”’*nd215-330cm* V"' s"‘, respec- 
tively. The resistivity of the seml-insuIating material de- 
creases gradually from about 1.5X10*/? cm toabout 5x10^ 
/? on as the As concentration in the melt increases from 
0.475 to about 0.535. The Hall mobility increases from about 
! 1500cm* V“' 8“‘ to about 6000 cm* V”’ s~',andthefree- 

electron concentration increases from about 2 X 10* cm~* to 
: about 2.5X10* cm~* over the same melt composition range. 

I The resistivity decreases sharply above about 0.535 As atom 
fraction as a result of an increase in the free-electron 
concentration. 

Detailed resistivity profiles of crystals grown under As- 
’ rich and Ga-rich conditions, as shown in Fig. 2, further em- 

: phasize the importance of mslt stoichiometry in controlling 

1 the electrical compensation. As-rich (Ga-rich) melts become 
I progressively more As-rich‘ (Ga-rich) as the crystal is pulled 
I from the melt. Crystals grown from As-rich melts are semi- 
I insulating (n-type) uniformly from front to tail as long as the 
i As atom fraction is less than about 0.535. On the other hand, 
i crystals grown under Ga-rich conditions below the critical 
I composition are invariably p type (Fig. 2). Furthermore, 
crystals grown from Ga-rich melts initially above the critical 
I composition undergo a transition from semi-insulating to p 
i type when As in the melt becomes depleted and the composi- 
\ tion reaches the critical value. This behavior clearly shows 
• that the semi-insulating to p-type transition is controlled by 
the melt stoichiometry and does not occur as a result of the 
! normal segregation of background impurities in the melt. In 
' addition, our analysis of the distribution of carbon along the 
crystals by LVM showed no correlation between carbon seg- 
regation and stoichiometry. 

An absorption band between 1 and 1.4 pm previously 
identified* with the EL2 center was observed in all of our 
semi-insulating material. The activation energy of the deep 
level as determined from variable-temperature Hall mea- 
surements* ( ~0.76 eV), and the behavior of the photocon- 



FIC. 2. Resistivity profiles of crysuk grown under As- and G«-rich condi- 
• tions. The initisl melt compositions for the As-rich, Gs-rich (As initially 
above critical composition), and Ga-rich (As initially below critical compo- 
sition) crystals are 0.501, 0,481, and 0.462 As atom fraction, respectively. 
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FlGs 3, Concentration of EL2 as deteimined by optical absorption mea- 
surements as a function of melt composition » llie concentration of EL2 
increases from 5x10** cm"^ to 1.7 X 10** cm“* as the As atom fraction 
increases from 0.48 to 0.5 1 and appean to saturate as the As concentration 
increases further to 0.53. 


ductivity* thresholds above and below 120 K were also 
found to be consistent with the presence of EL2. The concen- 
tration of EL2 was determined from the absorption using the 
cross section reported by Martin et a/.’ Absorption due to 
unoccupied EL2 centers was not observed, and the Hall 
measurements indicated that the centers were always more 
than 90% occupied. Consequently, the absoprtion was taken 
to be proportional to the total concentration of EL2. The 
concentration of EL2 was found to depend on the melt stoi- 
chiometry, as shown in Fig. 3, increasing from about 5 X 10'* 
cm~* to 1.7 X 10'^ cm~* as the As atom fraction increased 
from 0.48 to 0.51. The concentration remained constant as 
the As fraction increased further to about 0.53. These results 



FIG. 4. Carbon concentration as a function of the ratio of the EL2 concen- 
tration to the electron concentration evaluated according to the expression 
given in the legend. The linear behavior indicates the dominant role of EL2 
and carbon in controlling the electrical compensation. 
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•re consistent with data obtained on Retype material grown 
by vapor phase epitaxy and organometallic chemical vapor 
deposition." 

Using steady-state generation-recombination consider- 
ations we can obtain the expression — No 
(e7-/<7V7-), where iV^ and A'’/) are the concentrations of accep- 
tors and donors (other than EL2), is the concentration of 
occupied deep donors, n is the free-carrier concentration, e-f 
(Ref. 12) is the thermal emission rate, a (Ref. 13) is the cap- 
ture cross section, and Vr is the thermal velocity. We evalu- 
ated our results accordingly by plotting as a function of 
Ny/n taking and n as the concentration of carbon, 
occupied EL2 centers, and free-electrons, respectively. The 
concentrations of carbon (by LVM), EL2 and free-electrons 
were individually measured for each sample. The plot is lin- 
ear, as shown in Fig. 4, clearly indicating the dominant role 
of EL2 and carbon in controlling the electrical compensa- 
tion. The fluctuations in these data could be attributed to 
variations in the concentration of other background impuri- 
ties. The slope, of the curve in Fig. 4 is about a factor of 10 
higher than the value calculated using published values of 
the cross section and emission rate for EL2. The discrepancy 
could be due to errors in the published values for the optical 
cross sections of carbon and EL2. Another possibility is that 
carbon behaves either as an amphoteric or inactive impurity 
in LEC material so that the activation of carbon acceptors is 
only a fraction of the total concentration. 

In summary, we have shown that the melt stoichio- 
metry controls the incorporation of the EL2 center in un- 
doped LEC GaAs. The EL2 center, in turn, controls the 
semi-insulating properties of the material by compensating 
residual acceptors, which are principally carbon, The results 
suggest that the identity of EL2 is associated with an intrin- 
sic (native) defect such as an arsenic interstitial or a vacant 
gallium site. However, we cannot rule out the possibility that 
the segregation of oxygen, an impurity that has been associ- 
ated with high resistivity in GaAs,"''* depends strongly on 


the stoichiometry. In practical terms, our results show that 
an essential element in achieving uniform, high resistivity 
crystals is the close control of the As atom fraction in the 
melt. A more detailed account of this work will be presented 
in a furture paper. 

The authors are grateful to Dr. Peter A.sbeck for many 
helpful discussions during the course of this work. We also 
appreciate the assistance of Dr. Jane Yang in making the 
Hall measurements. We also wish to thank the National 
Aeronautics and Space Administration for partially sup- 
porting (his work under ccmtract No. NAS3-22224. 
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ABSTRACT I 


It Is shown that the key to reproducible growth of undoped semi- 
Insulatlng GaAs by the liquid encapsulated Czochralski technique 1s the control 
over the melt stoichiometry. Twelve crystals were grown from stoichiometric and 
non-stolchiometric melts. The material was characterized by secondary ion mass 
spectrometry, localized vibrational mode far Infrared spectroscopy. Hall effect 
measurements, optical absorption, and photoluminescence. A quantitative model 
for the compensation mechanism in the semi -insulating material was developed 
based on these measurements. The free carrier concentration is controlled by the 
balance between EL2 deep donors and carbon acceptors; furthermore, the incorpora- 
tion of EL2 Is controlled by the melt stoichiometry, Increasing as the As atom 
fraction in the melt Increases. As a result, semi -insulating material can be 
grown only from melts above a critical As composition. The practical signifi- 
cance of these results Is discussed in terms of achieving high yield and 
reproducibility In the crystal growth process. 
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A major Impediment to the development of GaAs Integrated circuits and 
microwave devices has been the lack of a reliable supply of semi -Insulating sub- 
strates with reproducible and thermal ly-stable [1] properties suitable for device 
fabrication. The liquid encapsulated Czochralskl (LEC) technique Is currently 
receiving considerable attention for the production of GaAs crystals because 
semi -insulating material can be grown without intentional [2] doping, and the 
technique offers the potential [3,4] for producing round, large-area substrates 
with uniform properties at a reasonable cost. 

The main question surrounding semi -insulating LEC GaAs has concerned the 
compensation mechanism by which the undoped material is semi -insulating. The 
understanding of the compensation mechanism has two important practical consequen- 
ces. First, knowledge of the cause-effect relationships between crystal growth and 
electrical characteristics of the material would greatly improve the yield of semi- 
insulating cystals in the growth process, as well as the crystal -to-crystal and 
wafer -to -wafer reproducibility. Second, this understanding would improve device 
performance. For example, backgating effects could possibly be diminished by 
adjusting [5] trap levels in material Intended for integrated circuit processing. 

It has been shown that the melt stoichiometry [6] is the key parameter 
which must be controlled to grow undoped semi -insulating LEC GaAs. This finding has 
led to the development of a quantitative model for the compensation mechanism. The 
purpose of this paper is to give a detailed account of these experimental results 
and interpretations. In the first section of the paper it is shown that the 
resistivity Is controlled by the melt stoichiometry. The second section describes 
the Identification of the predominant chemical impurities and electrically 
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The compensation model Is developed 


In the third section. 


II. APPAP.ATUS AND PROCEDURES 


The crystals used In this study were grown In the <100> direction from 
pyrolytic boron nitride (PBN) crucibles using a "Melbourn" high pressure (Metals 
Research, Ltd.) crystal puller. The pressure of the argon ambient during growth was 
300 psi. All of the crystals were undoped. The crystal weight and diameter were 
typically 2.5 Kg and 3 in., respectively. Nominally dry (< 500 ppm HgO) B 2 O 3 
encapsulating material (Johnson Matthey) was used for all growth experiments. The 
melt was prepared by in situ [7] synthesis, to achieve high purity, starting with a 
charge of 6-9' s As (Cominco) and 6-9' s 6 a (Kawecki). The composition of the melt 
was varied by changing the relative amounts of Ga and As in the initial charge. 

The composition of the initial melt was determined by adjusting the charge 
composition for the amount of As lost [ 6 ] during the heat-up cycle. The weight 
of lost As was determined by comparing the weight of the initial charge to the 
weight of the crystal and the charge remaining in the crucible after growth. The 
B 2 O 3 was carefully separated from the charge remaining in the crucible before 
weighing to eliminate possible errors introduced by B 2 O 3 evaporation. Samples 
for characterization were obtained from along the length of each crystal. The 
melt composition corresponding to the growth of each sample was determined by 
adjusting the initial melt composition for the crystal weight at the time of 
growth. The crystal weight and length were recorded during the growth as a 
function of time. After terminating the growth process, the crystals were cooled 
at between 30 and 80°C hr"^. 

The concentration of background transition metal and group IIl-VI impur- 
ities was determined by secondary ion mass spectrometry [ 8 ] (SIMS). The 
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concentration of carbon (acceptors) Mas determined by localized vibrational mode 
(LVM) far-infrared [9] spectroscopy. The LVM absorption band for occurs at 
582 cm"^. The Integrated absorption Is proportional to the concentration. 

The carbon concentration of our material was determined by measuring the 
absorption at 77K and converting to concentration by using the optical cross- 
section given in Ref, [9], The free carrier concentration. Hall mobility, and 
resistivity were determined by Van der Pauw-Hall effect measurements at room 
temperature and over a temperature range of from 298 to 400K, Optical absorption 
measurements In the near infrared were made at room temperature with a Perkln- 
Elmer Model 330 spectrometer, Photo! umlnescence measurements [10] were made at 
4,2K by Immersing samples In liquid helium, 

ni. STOICHIOMETRY-CONTROLLED RESISTIVITY 

The electrical characteristics of the crystals were evaluated by Hall- 
effect measurements using samples obtained from the fronts and tails of the 
twelve crystals, and from detailed resistivity profiles of five crystals. The 
resistivity was found to be a strong function of the melt stoichiometry, as shown 
in Fig. 1. 

Figure 1 shows that the material is semt-lnsulati ng (n-type) above, and 
p-type (low resistivity) below a critical As concentration in the melt of about 
0.475 As atom fraction. The resistivity peaks at the critical composition at a 
value of about 1.5 x 10® n-cm and decreases approximately 8 orders of magnitude 
below the critical composition. The resistivity also decreases very gradually as 
the As fraction increases from the critical composition. 

The variation in resistivity across the melt composition range is 
explained in terms of the corresponding free carrier concentration and Hall 
mobility, as shown in Figs, 2 and 3, The semi -i nsul ati ng material grown at the 
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critical composition is n-type with a carrier concentration and mobility of 
1-2 X 10^ cm’^ and 1-2 x lo^ cm^ V'^ s*^, respectively. These Hall mobilities 
are low for n-type GaAs possibly indicating the presence of inhomogeneities in 
the material. As the As atom fraction increases from the critical composition to 
about 0.51, the mobility gradually Increases to 4-5 x 10^ cm^ s*^, which is 
more typical of n-type material. The corresponding electron concentration 
gradually increases to 6-8 x 10^ cm’^. The combined increase of both the 
mobility and carrier concentration lead to a reduction in resistivity of about 
one order of magnitude, A further change in the As fraction from 0.51 to about 
0.54 brings about a more pronounced increase in free carrier concentration of 
about 2 orders of magnitude. The result is a corresponding decrease in 
resistivity to as low as 1.5 x lo® P-cm. 

The material becomes p-type below the critical composition. The free 
hole concentration rises approximately 9 orders of magnitude following a IS 
reduction in As fraction in the melt from the critical composition. The hole 
concentration and Hall mobility of this material are in the range of 
1-3 X 10^® cm"3 and 215-330 cm^ V"^ s"^, respectively. Mobilities obtained from 
the p-type material grown in the transition region, corresponding to melt compo- 
sitions within about 1% of the critical composition, were very low, between 1 and 
30 cm^ s”^. The measured hole concentrations were between about 1 x lo*^ and 
2 X 10l4 cm"^. These carrier concentrations are too high to explain the low 
mobilities in terms of mixed conduction. The low mobilities of material grown in 
the transition region could reflect inhomogeneities in the material. For 
instance, a striated pattern of regions of high and low resistivity could cause 
such behavior. 

Detailed resistivity profiles of crystals grown from initially As- and 
Ga-rich melts further emphasize the role of the melt stoichiometry in controlling 
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the electrical compensation. It is Important to note that unless the Initial 
melt Is precisely stoichiometric [11], As-rich (Ga-rich) melts become progres- 
sively more As-rich (Ga-rich) as the crystal is pulled from the melt. Crystals 
grown from As-rich melts were Invariably seml-lnsulatlng from front to tall. 
Crystals grown from 6a-r1ch melts Initially below the critical composition were 
p-type throughout. On the other hand, crystals grown from Ga-rich melts 
Initially above the critical composition underwent a transition from semi- 
insulating to p-type at the point along the crystal where the corresponding melt 
composition reached the critical value. This behavior clearly indicates that the 
resistivity is controlled by the melt stoichiometry, and, that the semi- 
insulatlng-to-p-type transition is not related to the normal segregation of some 
common background Impurity toward the tall of the crystal. Otherwise, the tall 
of As -rich-grown crystals would have become p-type as well. 


IV. ANALYSIS OF CHEMICAL IMPURITIES AND ELECTRICALLY ACTIVE CENTERS 

To establish a model for the compensation mechanism In undoped semi- 
Insulating LEC GaAs, the chemical impurities and electrically active centers were 
evaluated. The chemical impurities, as determined by SIMS and LVM, are shown in 
Table I. These results represent an average of measurements obtained from the 
front and tail of each of the twelve crystals. A range of concentrations is 
reported for those elements which exhibited substantial concentration variations 
from crystal to crystal. Carbon, as determined by LVM, Is the predominant, 
chemically identifiable, electrically active background impurity. The concentra- 
tion of carbon ranged from 2 * 10^^ to 1.3 x cm'^. The concentration of the 
ofher background impurities common to GaAs, such as S, Mg* Cr, Mn, and Fe, Is 
very low. In particular, the Si concentration Is consistently less than about 
1 X 10^® cm‘^ as determined by SIMS measurements. Low Si contamination results 
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from crystal growth without a quartz crucible. Si levels In this material are 
more than one order of magnitude lower than In Bridgman material that was 
analyzed. Boron Is present In a range from about 1 * 10^® to 6 x 10^^ cm”^. The 
source of boron contamination Is the B 2 O 3 encapsulant. The Incorporation of 
boron has been shown £ 12 ] to depend on the water content of the encapsulant. 
decreasing as the water content Increases. Although boron Is the predominant 
chemical Impurity. It Is Iso-electronic with Ga. and no evidence was found In 
this study Indicating that boron Is electrically active. 

The evaluation of the electrical and optical properties of the seml- 
insulatlng material Indicates that the deep donor, commonly referred to as EL2. 

Is the predominant deep center. An optical absorption band between 1 and 
1.4 urn previously identified with the EL2 center £13] was observed In all of the 
semi -Insulating material. In addition, the activation energy of the electron 
concentration, obtained from plots of the temperature-corrected free electron 
concentration as a function of the reciprocal of temperature, was 0.75 + 0.02 eV. 

This energy is consistent with published values [14] for the activation energy of 
EL 2 . The behavior of the photoconductivity thresholds [15] above and below 120 K 
was also found to be consistent with the presence of EL2. 

The concentration of EL2 was determined by optical absorption using the 
cross-section reported by Martin et al . [16]. Absorption due to unoccupied EL2 
centers was not observed, and the Hall measurements indicated that the center was 
more than 90% occupied. Consequently, the absorption was taken to be proportion- 
al to the total EL2 concentration. The concentration of EL2 was found to depend 
on the melt stoichiometry, as shown in Fig. 4. increasing from about 5 x 10^® 
to 1.7 X 10^^ cm”^ as the As atom fraction increased from about 0.48 to 0.51. 

The concentration remained constant as the As fraction Increased further to about 
0.535. 


49 


ORiQiNAL PAGE IS 
OF POOR QUALITY 



R<‘>ckwell International 


Photolunilnescence (PL) studies of the semi •insulating material are 
consistent with the measured dependence of EL2 on melt stoichiometry. Typical PL 
spectra of semi -insulating material grown from Ga- and As-rich melts, shown in 
Fig. 5, exhibit bands peaking at 0.68 and 0.77 eV. The 0.68 eV band has been 
attributed [10] to radiative-recombination between EL2 electron traps and the 
valence band, and the 0.77 eV band to recombination possibly associated [10] with 
a hole trap. The intensity of the 0.68 eV band in the semi -insulating GaA'^^ grown 
from Ga-rich melts is substantially reduced by comparison with As-rich grown 
material. This behavior is consistent with the decrease of the EL2 concentration 
with decreasing As fraction (Fig. 4) as determined by optical absorption. 

Neither band was observed in the p-type material. 


V. THE COMPENSATION MECHANISM 

To develop a model for the electrical compensation in terms of the 
concentration of a predominant electrically active centers in the semi -insulating 
material, the concentration of shallow and deep centers were related through 
first-principle theoretical considerations. The ionization of EL2 produces an 
ionized center plus an electron in the conduction band; 

Unionized EL2 ^ Ionized EL2 + e' . (1) 


According to the principle of detailed balance, the concentration of ionized 
centers Nj, the concentration of electrons n, and the concentration of unionized 
centers Ny, are related by the following equation; 

N.n 

K , (2) 

where K is a constant determined by the thermodynamics of the system. Nj is 
equal to the net acceptor concentration, given as the difference in concentration 
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Nj - Nft - % , (3) 

The concentration of acceptors Is given as the surn of the concentrations of 

D 

carbon and other residual acceptors, N^; 

■ [carbon] + . (4) 

The concentration of unionized centers Is equal to the EL2 concentration as 
determined by optical absorption. That Is, only EL2 centers which are occupied 
by electrons contribute to the optical absorption process; 

Ny » EL2 . (5) 

By substituting Eqs. (3) through (5) Into Eq. (2), the following expression for 
the free-electron concentration Is obtained In terms of the predominant centers 
in the material; 


n 


EL2 

([carbon] + nJ - Ny) 


This expression can be rewritten In the following form: 
[carbon] * K 


(6) 


(7) 


Therefore, the carbon concentration Is proportional to the ratio of the EL2 
concentration to the electron concentration. 

Semi-Insulating material grown from melts ranging from 0.475 to 0.535 As 
atom fraction was evaluated according to Eq. (7) from measurements of the carbon 
concentration (by LVM), the EL2 concentration (by optical absorption), and the 
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electron concentration (by Hall effect measurements) for each sample. A plot of 
the carbon concentration as a function of the ratio of the EL2 concentration to 
the electron concentration, shown in Fig. 6, follows linear behavior, indicating 
that the electron concentration is indeed controlled by the balance between EL2 
and carbon. This result is independent of possible errors in the '•ublished 
values of the optical cross-sections for carbon and EL2. It is important to note 
that if some other impurity were the predominant acceptor, such as Mn, Fe, Cu, or 
Zn, the linearity predicted on the basis of Eq. (7) would still necessarily 
hold. However, the linearity would not be distinguishable because the term 
[carbon] would be small compared to N^; the figure would -be a scatter plot. In 
fact, the scatter in these data probably reflects actual fluctuations in the 
concentration of other background impurities rather than random error in the 
experimental measurements. The small value of the intercept (Nq - N^) of the 
least-squares fit to the data also indicates the predominance of carbon accep- 
tors. Thus, EL2 deep donors and carbon acceptors control the elecrical compen- 
sation in semi -insulating LEC GaAs grown from melts ranging from 0.475 to 0.535 
As atom fraction. 


The variation of the electrical characteristics of the semi -insul ating 
material (Figs. 1-3) with melt stoichiometry can now be explained on the basis of 
the analysis above. The EL2 concentration must either exactly match or exceed 
the carbon concentration to produce semi -insulating properties. The EL2 concen- 
tration in material grown from Ga-rich melts below the critical composition is 
insufficient to compensate the carbon, leading to p-type conductivity. Semi- 
insulating material grown at the critical composition is closely compensated, 
leading to the maxima of the resistivity. As the As atom fraction in the melt 
increases from the critical composition to about 0.51, the EL2 concentration 
becomes progressively higher than the carbon concentation. As a result, thermal 
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ionization of (unionized) EL2 centers (see Eq. (1)) gives rise to a gradual 
increase in the electron concentration and a corresponding decrease in the 
resistivity (Figs. I and 2). Our model further shows that the relatively low 
resistivity of semi -insulating material grown from melts between 0,53 and 0.54 As 
atom fraction is due to exceptionally lew concentrations of carbon acceptors. 


VI. SUMMARY AND CONCLUSIONS 


It was shown that the melt stoichiometry controls the compensation in 
undoped 6aAs grown by the liquid encapsulated Czochral ski technique. The free 
carrier concentration of the semi -insulating material is controlled by the 
balance between EL2 donors and carbon acceptors. The incorporation of E12 
depends on the melt stoichiometry, increasing as the As melt concentration 
increases. As a result, semi -insulating material can be grown only above a 
critical As concentration. 


From a fundamental standpoint, these results have provided the first 
direct evidence indicating that the identity of the EL2 center is associated with 
a native defect or native-defect complex in bulk GaAs: LEG QaAs is pure enough 

to reveal the contribution of native-defects in bulk material. The dependence of 
EL2 on melt stoichiometry is consistent with published work on GaAs grown by 
vapor phase epitaxy [173 and organometallic chemical vapor deposition [18]. 
Isolated native defects which would follow the stoichiometry dependence of ELE 
(Fig. 4) include the gallium vacancy Vqj,, the arsenic interstitial As^ , and the 
antisite Asq^. Since Vgj, would be expected to be an acceptor, ELE would more 
likely be related to one of the latter two defects. There is evidence [10,19-24] 
supporting either defect, and we believe that further work on welV characterized 
material is needed to make a definitive assignment of EL2. 
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In practical terns, these results shoM that setnl -Insulating GaAs can be 
grown by the LIC technique reproducibly and with high yield provided that the 
melt is As-rich* This condition ensures that the melt will not become 6a-r1ch 
during the growth process. Indeed, nine crystals grown during the course of this 
investigation from near-stoichiometric As-rich nelts were semi -Insulating from 
front to tail. In addition, the characteristics of Se active layers Implanted 
directly into the substrates showed excellent uniformity [25] of the pinch-off 
voltage across each wafer. These results demonstrate that high yield and repro- 
ducibility can be achieved In both the growth of GaAs by the LEG technique and in 
the application of LEG GaAs to the Integrated circuit technology. 
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Fig, 1 Dependence of electrical resistivity of LEG GaAs on melt stoichiometry. 

Semi -insulating material is obtained above, and p-type (loM-resistivity) 
material below a critical melt composition of about 0.475 As atom 
fraction. 

Fig. 2 Dependence of free carrier concentration of LEG GaAs on melt 

stoichiometry. The semi -insulating material is n-type, and the free 
electron concentration increases gradually as the As concentration in 
the melt increases from the critical composition. The material becomes 
p-type below the critical composition. The free carrier concentration 
rises approximately 9 orders of magnitude following a 1% reduction in As 
fraction from the critical composition. 

Fig. 3 Dependence of Hall mobility of LEG GaAs on melt stoichiometry. The 
mobility of the semi -insulating material varies from 1-2 x 10^ to 
4-5 X 10^ cm^ V"^ s"^ as the As atom fraction increases from the cri- 
tical composition to about 0,535. The mobility of the p-type material 
grown in the transition region, within about 1% of the critical composi- 
tion, is low, between 1-30 cm^ s"^. The mobility of the p-type 
material grown outside of the transition region ranges from 215 to 330 
cm2 s*^ 


Fig. 4 
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Dependence of EL2 concentration as determined by optical absorption on 
the melt stoichiometry. The concentration of EL2 Increases from 
5 X 10^® cm'^ to 1.7 X 10^® cm'3 as the As atom fraction Increases from 
about 0.48 to 0.51. and appears to saturate as the As concentration 
Increases further to 0,53. 


Fig. 5 Typical photoluminescence spectra of semi -Insulating GaAs grown from (a) 
an As-rich melt (As atom fraction ■ 0.507, p ■ 1.8 x lO^ n-cm) and (b) a 
Ga-rich melt (As atom fraction « 0.438, p » 1.4 x lo® o-cm). The 
Intensity of the 0.68 eV band decreases as the As atom fraction 
decreases toward the critical composition consistent with the optical 
absorption measurements (Fig. 4). 


Fig. 6 Dependence of the carbon concentration on the ratio of the EL2 con- 
centration to the electron concentration. The concentration of carbon, 
EL2, and electrons was detennined for each sample. The dashed line Is a 
least-square fit to the data. The linearity of the data Indicates the 
dominant roles played by EL2 deep donors and carbon acceptors In 
controlling the compensation (see text). The small value of the 
intercept (Nq - N^) also Indicates the predominance of carbon acceptor. 
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Table I. Chemical Impurity analysis of 12 LEC GaAs crystals. 
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DISLOCATION STUDIES IN 3-INCH DIAMETER 
. LIQUID ENCAPSULATED CZOCKRALSKI GaAs 

R. T. Chen and D. E. Holmes 
Rockwell International 

Microelectronics Research and Development Center 
Thousand Oaks* California 91360 

I i ■ 

ABSTRACT f 

He have characterized the density and distribution of dislocations In 
3-1nch diameter, undoped GaAs crystals grown by the liquid encapsulated 

i 

Czochralski technique. The radial distribution across wafers follows a 
"W'-shaped profile Indicating excessive thermal gradient- Induced stress as the 
primary cause of dislocations. The density along the body of each crystal 
Increases continuously from front to tail. In contrast, the longitudinal 
distribution In the cone region Is Inverted, first Increasing, and then 
decreasing as the crystal expands from the neck to full diameter. Growth 
parameters favoring reduced dislocation densities include good diameter 
control, the use of thick B 2 O 3 encapsulating layers, slightly As-rich melts, 
and low ambient pressures. The dislocation density in the body of the crystal 
Is practically Independent of cone angle e for 20* < e < 70*. However, high 
densities result for flat-top (0* < e < 20*) crystals. Dash-type seed necking 
works to reduce the dislocation density only when high-density seeds 
(>5000 cm“^) are used. Dislocation densities below 1 x 10^ cw"^ can be 
achieved routinely In the front of a crystal through proper control of these 
parameters. The lowest measured dislocation density was 6000 cm“^. Further, 
we show that convective heat transfer from the crystal to the high pressure 
ambient plays a dominant role in controlling the dislocation density. 
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1. INTRODUCTION 

Current interest in large-diameter GaAs crystals grown by the liquid 
encapsulated Czochralski (LEC) technique [1,2] stems from the need for uni- 
formly round substrates with stable, semi-insulating properties for digital 
and monolithic integrated circuits. As these circuits become larger and more 
complex, adverse effects of dislocations on device performance and reliability 
could appear. There have been virtually no systematic studies reported con- 
cerning the possible role of dislocations because of the early stage of devel- 
opment of GaAs integrated circuits. It is therefore important to characterize 
the density and distribution of dislocations across large-diameter substrates, 
and to establish how the density and distribution is controlled by the crystal 
growth process. This understanding is also important for the application of 
large-diameter LEG material to minority-carrier devices such as solar cells, 
where low dislocation densities are required to achieve high minority-carrier 
lifetimes and diffusion lengths. 

A principal cause of dislocations in bulk GaAs crystals is stress 
induced by thermal gradients [3-6] during crystal growth. Radial gradients 
are of particular concern in Czochralski -type growth configurations (LEG, 
Grenmelmaire). Most of the published dislocation studies on GaAs concern 
small -diameter (<0.5 inch- diameter) crystals grown by the LEG [7-9], Bridgman 
[10,11], and modified Gremmelmaire [12] techniques. Since gradients generally 
decrease as the crystal diameter decreases, “effectively*' dislocation- free, 
small -diameter GaAs crystals have been grown [7,9,11,12]. Growth parameters 
reported to reduce radial gradients in small-diameter LEG crystals include the 
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height of the B2O3 encapsulating layer [7,13] and the cone angle [14]. 

Material properties which have been identified with the suppression of 
dislocations include the concentration of impurities [9,15.16] and melt 
stoichiometry [8,10,11]. To our knowledge, extensive dislocation studies of 
large-diameter LEG GaAs have not been published. 

i 

In this paper, we present and discuss results of dislocation studies 
involving undoped, 3- inch- diameter LEG GaAs crystals. We investigated the 
radial distribution across substrates, as well as the longitudinal distribu- 
tion along the length of the crystals. The effects of cone angle, B2O3 encap- 
sulant thickness, ambient pressure, seed quality. Dash- type necking, diameter 
control, and melt stoichiometry were evaluated quantitatively. As a result of 
this effort, we showed that dislocation densities below 10^ cm"^ can be 
achieved routinely through proper control of the crystal growth process. 
Dislocation densities as low as 6000 cm"^ are reported. Finally, techniques 
for possible further reductions in dislocation density are discussed. 


II. APPARATUS AND PROGEDURES 

Eighteen undoped, 3-inch diameter GaAs crystals were grown during this 
study using a high pressure "Melbourn" (Metals Research Ltd.) puller. The 
crystals were grown in the <100> direction from 6-inch-diameter quartz and 
pyrolytic boron nitride (PBN) crucibles. The B2O3 encapsulant was obtained 
from Johnson-Mat they. The seeds were cut from Bridgman crystals and had 
dimensions ofAimxAmmx 50 mm. The GaAs melt was prepared by in situ 
synthesis [17] starting with a charge of high purity Ga (Kawecki) and As 
(Gominco) weighing 3 kg. The crystals weighed from 2.2 to 2.6 kg. The seed 
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•nd crucible were both rotated counter-clockwise at 6 and IS rpm» 
respectively. The crystal pull rate was 7 mm hr*^, and the crucible lift rate 
was 1.4 nm hr'^. When the growth process was terminated, the crystal was 
positioned above the B 2 O 3 encapsulating layer, and the system was cooled at a 
constant rate of between 30 and 80 *C hr”^. 

The diameter of the crystals was controlled manually by adjusting the 
temperature and cooling rate. Changes In temperature or cooling rate were 
made In response to variations of the differential weight of the crystal. The 
differential weight signal was provided by a "load cell," a weighing device on 
which the pull shaft was mounted. An Increase or decrease of the differential 
weight Indicated a corresponding expansion or contraction of the crystal. 

Growth was typically Initiated by reducing the diameter of the crystal below 
that of the seed. This procedure Is referred to as Dash- type seed necking 
[18]. The diameter was then allowed to expand controllably, forming the 
"cone." Expansion was terminated at the "shoulder," and the diameter of the 
crystal was held constant for the remainder of the growth process, forming the 
"body". Crucible lift was Introduced during the growth of the shoulder. 

The crystals were cut according to the diagram shown In Fig. 1. The 
samples were lapped and polished on both sides. Dislocation densities and 
distributions were evaluated by preferential etching (KOH for 25 min. at 400"C). 
This etch attacks dislocations Intersecting the surface of the sample, forming 
hexagonal etch pits. The etch pit density (EPO) corresponds directly to the 
dislocation density as confirmed by x-ray topography at our laboratory and 
elsewhere [19]. The EPO measurements were typically made from low-magnifica- 
tion (70X) micrographs by counting the pits over 1.3 x 1.0-m1111meter 
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regions. However, higher magnification (either 140X or 280X) was required to 
resolve the pits when the EPD exceeded approximately 1 x 10® cm"^. The esti- 
mated error In counting the pits on each micrograph was less than 5%. 

We investigated the effect of seven growth parameters on the disloca- 
tion density and distribution. These were: (1) cone angle, (2) seed quality, 
(3) seed necking, (4) diameter control, (5) melt stoichiometry, (6) height of 
the B2O3 encapsulating layer, and (7) ambient pressure. The cone angle, 
defined as the angle between the wall of the cone and the horizontal, was 
varied from 0 to 65“. Crystals with a cone angle less than 20“ are referred 
to as "flat-top" crystals. Referring to seed quality, the EPD of the seeds 
ranged from about 1.5 x 10^ to 5 x 10® cm"^. Crystals were grown with high 
and low EPD seeds with and without Dash- type necking. The neck diameter 
varied from 1.2 to 3 mm. Diameter control refers to the deviation of the 
diameter from the average value. We achieved our lowest diameter deviations 
(il.l nm) by controlling the diameter with the cooling rate with minimal 
(direct) adjustments of the temperature. The Initial melt stoichiometry was 
varied from 0.462 to 0.506 As atom fraction. Procedures necessary for making 
an accurate determination [20] of the melt composition have already been 
discussed. The height of the B2O3 encapsulant above the melt was approx- 
imately 17 mn In the majority of the growth experiments, corresponding to 
500 g of material. One experiment each was made with 170 and 390 g of B2O3, 
corresponding to 9 and 13 nm heights, respectively. The ambient pressure 
during growth was typically 300 psi. One experiment was conducted at a lower 
pressure of 50 psi. 
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I 
i 

111, results and discussion 

3.1 Radial Dislocation Distribution 

The distribution of dislocations across wafers exhibited four-fold 
symmetry indicative of the <100> crystallographic orientation, as shown in 

Fig. 2. A microscopic view of the dislocation distribution across wafers, as t 

I 

shown in Fig. 3, clearly shows large variations in EPD. The main features of 

I 

the distribution are; (1) minimum EPO occurs over a large annulus between the I 

il ■ 

i 

center and edge (Region 1, or “ring" region); (2) intermed'.ate EPD occurs in I 

I 

the center (Region 2, or "center" region); (3) maximum EPD occurs at the edge 
(Region 3, or "edge" region). In addition, the EPD in the ring and edge 
regions is greater along the <i00> than the <U0> direction (see Fig, 3). 

Measured EPD distributions across the full diameter of wafers typically 
followed a "Vi"-shaped profile, as shown in Fig. 4. 

Our experimentally determined radial EPD distributions are consistent 
with theoretical thermoelastic analyses of Czochralski crystals of Penning [3] 
and Jordan et al. [5,6]. Jordan calculated the total stress in the crystal in 
terms of twelve (111) <110> slip systems. The dislocation density is assumed 
to be proportional to the total stress within an additive constant. Since the 
periphery of the crystal is cooler than the center as the crystal is pulled 
from the melt, the periphery and center are under tension and compression, 
respectively. The calculated stress is highest at the periphery, consistent 
with the experimental finding that the maximum EPD occurs in the "edge" region 
of our crystals. The calculated stress is lowest in the transition between 
regions of tension and compression, consistent with the fact that the lowest 
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nteasured EPDs occur in the “ring" region. The relatively high EPD measured 
along <100> compared to along <1I0> is explained in terms of the theory by the 
fact that more slip systems contribute to the total stress along <100>. The 
agreement between theory and experiment indicates that radial gradient- induced 
stress is the principal cause of dislocations in our crystals. Our results 
also agree with other experimental studies [5] of (100) GaAs. 

Variations in the nature of the radial distribution from the front to 
the tail of the crystal are discussed in Section 3.3. 

3.2 Dislocation Networks; Polygonization 

An important inorphological feature of the microscopic dislocation dis- 
tribution is the formation of etch pit networks. Two types of morphologies 
are observed. The first is a cellular network, shown in Figs. 2 and 3, where 
the dislocations form an interconnected network of cells with few dislocations 
within each cell. The approximate diameter of the cells is 500 ym, corres- 
ponding to an EPD of 2 X 10^ cm“^. The cell diameter decreases as the EPO 
increases (100 ym diameter corresponds to an EPD of about 1 x 10® cm"^). When 
the EPD is less than about 2 x 10^ cm“^, the morphology of the network takes 
on a lineage structure, where the etch pits form visible wavy lines. These 
lines extend from a few millimeters to more than one centimeter, and are 
oriented along <110>, as shown in Fig. 2. 

We believe that these dislocation networks form as a result of the 
polygonization process [21], where the dislocations re-align after solidifica- 
tion to reduce the strain energy of the crystal. The re-alignment probably 
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occurs by both climb and glide processes. In general, dislocations in zinc 
blend materials can undergo alignment into walls defined by (110) planes 
perpendicular to (111) slip planes. These walls would intersect (100) planes 
along <110> directions, consistent with our observations. The cellular net- 
work, in effect, constitutes a high packing density of lineage structures 
where one dislocation may interact with several dislocation lines, forming 
Interconnected networks, 

3.3 Longitudinal Dislocation Distribution 

The longitudinal variation (Along the growth direction) of the disloca- 
tion density was examined by comparing radial distributions of wafers obtained 
from the front, middle, and tail of the crystals, as shown in Fig. 4. The EPD 
invariably increased from front to tail in each of the three regions, as shown 
in Table I, while the radial profiles remained "W-shaped, as shown In 
Fig. 4. This behavior could indicate that the overall level of stress 
increased along the crystal, or that the dislocations multiplied after growth, 
or both. The average EPD increased from front to tail (see Table I) by a 
factor of 8, 7, and 1.5 in the "ring," "center," and "edge" regions, respec- 
tively. Further, the ratio of the EPD in the "center" to that in the "ring" 
region decreased from front to tail in the majority of the crystals, as shown 
in the table. These results show that the radial EPD distribution becomes 
more uniform toward tail of the crystals even though the "W"-shaped pro- 
file persists. 
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3.4 Parameters Affecting Dislocation Density 

In this section, we present and discuss our results concerning the 
quantitative dependence of the dislocation density on cone angle, B 2 O 3 thick- 
ness. ambient pressure, seed quality and necking, diameter control, and melt 
stoichiometry. The effect of each growth parameter was evaluated by deter- 
mining the change of the EPD across each substrate as that parameter was 
independently varied. Since the EPO density in the “ring,** “center," and 
"edge" regions represented local limits of the entire EPO distribution, we 
were able to characterize the entire distribution with these three EPO 
values. Only when all three of these EPO values either increased or decreased 
were conclusions drawn concerning the effect of that particular growth 
parameter. We believe that the spatial resolution of our measurements In the 
center and ring regions (averaging over 1.3 x l.o im areas) was sufficient to 
reflect true variations in the average dislocation density across wafers while 
minimizing contributions due to microscopic fluctuations in density associated 
with polygonization. However, since higher magnifications were used to deter- 
mine the EPD near the edge of the crystals, we estimate that these measure- 
ments probably represent the true average EPD to within 50%. Therefore, 
measurements obtained from the center and ring regions were more sensitive 
indicators of actual EPD variations from crystal to crystal than measurements 
from the edge. In effect, we weighed the center and ring regions more 
heavily. (Since the center and ring measurements encompassed approximately 
80% of the area of a substrate, the heavier weighting of the center and ring 
measurements seems justifiable from a practical standpoint.) The EPD value 
reported in the tables is an average of at least two measurements. 


76 

J/3767A/cb 



Rockwell International 


3.4.1 Cone Angle 

We evaluated the effect of the cone angle on the dislocation density by 
comparing the EPD of full-diameter wafers cut from the front of each crystal. 
Our results, shown In Table II, show no correlation between cone angle and EPO 
for cone angles greater than about 25*. For example, crystals No. 9 and 
No. 10 were grown under very similar conditions In terms of the other six 
parameters evaluated in this paper. The only difference Is In the cone 
angles, which are 30* and 62* for crystals No. 10 and No. 9, respectively. 

The data shows virtually no difference between the EPD values In the center 
and ring regions. 

On the other hand, the EPD in the front of the flat- top crystal 
(No. 15) is In the low 10® cm"^ range, which is high compared to the other 
crystals. In addition, the longitudinal distribution is inverted along 
approximately the first half of the crystal, first decreasing from the front 
toward the tail before Increasing again as in all the other crystals. We also 
observed that the crystal began to expand rapidly when the top of the crystal 
emerged from the B 2 O 3 encapsulating layer, leaving a bulge at a distance from 
the front of the crystal equal to the height of the B 2 O 3 layer. This behavior 
shows that the crystal experienced significant additional cooling when 
emerging from the B 2 O 3 , indicating that the convective heat transfer from the 
crystal to the ambient was large compared to the heat transfer to the liquid 
encapsulating layer. The increased cooling presumably raised the level of 
stress near the top of the crystal leading to the unusually high dislocation 
density. 
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Dislocation maps of longitudinal cross sections of cones (see Fig. 5) 
were analyzed to follow the dislocation density distribution along the growth 
direction for various cone angles. The "VT-shaped radial dlstrlbutfon 
observed across wafers was clearly visible In these samples, as shown in- 
Fig. 5. However, the longitudinal EPD Increased after the neck, reached a 
maximum value, and then decreased before the crystal reached full diameter. 

(We had anticipated a continuous Increase In EPD In the cone region because 
the diameter expands continuously, and radial gradients typically Increase as 
the diameter Increases.) We found that the maximum value of the EPD decreases 
as the cone angle Increases, as Indicated In Fig. 5. A high concentration of 
slip traces was also observed In the cone region In crystals grown with low- 
angle cones. In addition, the maximum of the longitudinal EPD distribution 
was located directly below the neck In low-angle cones, and closer to the 
center of the high-angle cones, as Is evident In Fig. 5. 

The variation with cone angle of both the EPD at the maximum and the 
position of the maximum within the cone as the cone angle decreased from 65* 
to 30* Is consistent with the behavior of the flat-top crystal; i.e., the 
maximum EPD occurred at the top of the flat-top crystal, and the EPD at the 
maximum was the highest of all our crystals. We therefore believe that the 
same mechanism controlled the dislocation density and distribution at the top 
of all of the crystals, the flat- top crystal representing the limiting case of 
a 0* cone angle. In view of our discussion earlier In this section concerning 
the flat-top crystal, we further believe that the dislocation maximum formed 
as the cone emerged from the 6^03 encapsulating layer as a result of increased 
convective heat transfer to the ambient. The dislocations associated with the 
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maxlinuffl represent a ''secondary" distribution added to the primary ("grown- In") 
distribution, which formed at the solidification front. 

Our study of the secondary dislocation distribution suggests the 
following model for the heat flow In the crystal at a position corresponding 
to the top surface of the encapsulating layer. The Isotherm shape Is deter- 
mined by the relative vertical and radial components of heat flow. The verti- 
cal heat flow Is relatively large when the crystal Is thin (as the neck 
emerges from the B 2 O 3 ), and the Isotherm shape Is relatively flat. When the 
cone begins to emerge from the B 2 O 3 , radial heat flow becomes more Important; 
the Isotherm shape becomes more concave with respect to the solid as the 
radial gradient Increases. The radial gradient Increases as the cone angle 
decreases, leading to more pronounced EPD maximums for luw-angle cones. As 
the vertical wall of the crystal begins to emerge from the encapsulating 
layer, the curvature of the Isotherms and the radial gradients decrease. 

3.4.2 B2O3 Height 

We evaluated the effect of the height of the B 2 O 3 encapsulating layer 
on the EPD by varying the height from 9 to 17 mm. The results In Table III 
show that the EPD decreases as the height of the layer Increases. The effect 
Is more pronounced at the front of the crystals. In addition, the nature of 
the secondary dislocation distribution In the cone region was Independent of 
B 2 O 3 height. This behavior Indicates that the radial gradients near the 
crystal -melt Interface decrease as a direct result of the presence of a 
thicker B 2 O 3 layer. In view of the results of the previous section, which 
showed that the heat transfer from the crystal to the ambient (above the B 2 O 3 ) 
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Is greater than the heat transfer to the B 2 O 3 liquid, we believe that the 
reduction of the radial gradient In the crystal attributed to thicker B 2 O 3 
layers results from more effective thermal Isolation between the region of the 
crystal near the melt Interface and the Ar ambient. This finding disagrees 
with theory [5], which predicts that the radial gradient would decrease as the 
B 2 O 3 height decreases. 

3.4.3 Ambient Pressure 

One crystal (No. 14) was grown at low pressure (50 psi). EPD measure- 
ments obtained from the front of the crystal are shown In Table IV. Excessive 
thermal degradation took place at the surface of the crystal due to the low 
ambient pressure. As a result, Ga droplets, which formed at the cone, 
thermally migrated through the crystal to the tall. The presence of the Qa in 
the crystal prevented the measurement of the EPD In the tall. The degrade^ 
tion, and subsequent loss of As from the cyrstal during growth, also prevented 
making an accurate determination of the melt stoichiometry. However, the 
electrical characteristics of the material Indicated that both the Initial and 
final melt compositions were within the As-rich range [20] similar to crystal 
No. 16. A comparison of the EPDs of the cyrstal grown at low pressure 
(No. 14) and crystal No. 16 shows that the EPD of crystal No. 14 was lower, as 
shown In Table IV, Indicating that the use of lower ambient pressures Is 
effective in reducing the EPD. In fact, the EPD of 6000 cm”^ In the ring 
region was the lowest value observed In this study. 

The results reported In Sections 3.4.1 and 3.4.2 of this paper Indicate 
the Importance of convective heat transfer via the ambient In controlling the 
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dislocation density. The heat transfer coefficient of the crystal -ambient 
surface is expected to increase as the square root of the pressure according 
to Jordan [5,6]. Therefore, we could have expected a reduction In heat 
transfer coefficient by no more than a factor of 2.5 by reducing the pressure 
from 300 to 50 psi. Our experimental finding of a 50% reduction in EPD Is 
consistent with the theoretical prediction. 

3.4.4 Seed Quality and Necking 

We determined the effectiveness of the seed quality and the Dash- type 
necking procedure In reducing the EPD by growing crystals from high-and-low 
EPD seeds with and without thin necks. The crystals were evaluated by com- 
paring the EPDs in the front of each crystal at full diameter. The results, 
given in Table V, show that low-EPD crystals (EPD <2.5 x 10^ cm"^) can be 
grown by employing low-EPD seeds with and without necking as well as by 
employing high-EPD seeds with necking. 

To understand the effect of seed necking, we studied longitudinal cross 
sections of crystals in the neck region. We could not directly observe the 
grown-in EPD in this region with neck diameters of less than about 2.5 mm 
because the neck region apparently deformed under the weight of the crystal, 
as shown in Fig. 6(a). However, dramatic reductions in EPD were observed for 
necks between about 2.5 and 3.5 inn In diameter, as shown in Fig. 6(b). These 
results indicate that the Dash- type necking procedure indeed works to reduce 
the dislocation density Independent of the EPD of the seeds. Yet, the effect 
was registered in the first full -diameter wafer only for high EPD seeds. We 
interpreted this behavior to mean that dislocations can be transmitted from 
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the seed to the crystal, and the transmission Is reduced by necking* However, 
the effect of necking Is limited; dislocations will be generated In the 
crystal even If the seed were perfectly dislocation free. 

3.4.5 Diameter Control 

The effect of diameter control was evaluated by comparing the EPD of 
crystal tails. Our results, shown In Table VI, Indicate that Improved 
diameter control generally reduces the dislocation density. For example, 
Crystals No. 6 and No. 9 were grown under very similar conditions, except that 
the diameter deviation was smaller In No. 6. The EPOs In the front of No. 6 
are higher than In the front of No. 9, whereas the EPDs In the tall are 
lower. We attributed the lower EPD In the tail of No. 6 to the Improved 
diameter control. We note that the effect of diameter control Is much less 
pronounced compared to that of cone angle, seed quality, and seed necking. We 
believe that crystals with more unstable diameter control were subjected to 
greater transient gradient-induced stress, which resulted In higher EPDs. 

3.4.6 Melt Stoichiometry 

We studied the effect of melt stoichiometry on the dislocation density 
by growing crystals from stoichiometric and non- stoichiometric melts. We 
found no correlation between EPD and melt stoichiometry for Ga- or As-rich 
melts with compositions less than 0.503 As atom fraction, as shown In 
Table VI. However, we note that the growth conditions and physical parameters 
of crystals No. 11 and No. 12 are nearly Identical except for the melt com- 
position. Yet, the EPD values In the front of crystal No. 11 are slgnlfl- 
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cantly lower compared to crystal No. 12. The reduced EPD values In the front 
of the crystal would Indicate that the As-rich melt favors reduced dislocation 
densities for melt compositions greater than aboutO.SOS As atom fraction. No 
significant Improvement Is apparent In the tall of No. 11, suggesting that a 
small range of melt compositions between 0.50S and 0.535 provides for optimal 
EPD reductions. 

IV. SUMMARY AND CONCLUSIONS 

We have characterized the density and distribution of dislocations In 
3-1 nch diameter LEC GaAs crystals. The radial distribution across wafers Is 
“W"-shaped, Indicating excessive thermal gradient- Induced stress as the pri- 
mary cause of dislocations as predicted on the basis of the models of Penning 
[3i] and Jordan et al . [5,6]. The dislocation density along the crystals 
Increases from front to tall at full diameter. Indicating that the level of 
stress In crystals grown In our current configuration Increases as the crystal 
Is pulled from the melt, or that the dislocations multiply after growth. The 
radial EPD distribution becomes more uniform toward the tall of the crystals 
even though the “W"-shaped distribution prevails. Jordan et al. [5] noted 
that a more "diffuse" radial dislocation distribution could result from the 
movement of 60-degree dislocations out of their slip planes Into the next-to- 
grpw layer of the crystal, adding to the glide dislocations at the solidifi- 
cation front. This explanation would seem to be a reasonable basis for 
modelling the observed behavior. However, such 'a model would have to take 
Into account the fact that the Increased uniformity of the dislocation density 
across 3-1nch substrates Is a very long-range process. 
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We have also evaluated the dependence of the dislocation density on 
seven crystal growth parameters. The EPD at full-diameter crystal Is 
virtually Independent of the cone angle e for 20* < e < 70*. However, the EPO 
Increases significantly for 0* < e < 20*. Our analaysis of the longitudinal 
dislocation distribution within the cone region further shows that the EPD 
Inverts for 20* < 0 < 70*, first Increasing, and then decreasing from front to 
tall before the crystal expands to full diameter. We explained these findings 
In terms of the dominant role played by convective heat transfer from the 
crystal to the ambient gas as the crystal emerges from the B 2 O 3 . In practical 
terms, these results show that crystals can be grown with a minimum disloca- 
tion density using a cone angle of about 30*. The use of a 30* cone angle 
maximizes the number of low-dislocation wafers that can be obtained from 
crystals while minimizing the time required to grow the cone. 

The dislocation density In tfie front of the crystals was found to be a 
relatively strong function of the height of the encapsulating layer, decreas- 
ing as the layer height Increases. We showed that this effect is probably a 
direct result of reduced radial gradients in the crystal within the B 2 O 3 layer 
due to more effective Insulation from the ambient gas. 

We showed that the Dash- type seed necking procedure is effective In 
reducing the dislocation density only when the EPD of the seed Is high 
(> 5000 cm"^); low dislocation crystals were grown with poor quality seeds 
with necking, and with high quality seeds with and without necking. These 
results Indicate that dislocations Indeed transmit from the seed to the 
crystal, and necking greatly reduces this effect. However, the mechanism can 
have only a limited effect, as dislocations are generated In the crystal even 
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if the seeds were perfectly dislocation- free. We would expect that further 
reductions, in the dislocation density in our crystals (achieved by altering 
the present thermal configuration) will require higher quality seeds if tiie 
necking procedure were to be eliminated from the growth process. The elim- 
ination of necking through careful selection of seeds would be advantageous 
from the practical standpoint of minimizing the time required to grow a 
crystal . 

We have also found that good diameter control and the use of slightly 
As-rich melts favor reduced dislocation densities. However, these effects are 
small compared to those of the cone angle, B2O3 height, seed quality and 
necking. Further work is needed to understand the effect of the melt 
stoichiometry. 

Finally, the experimental results presented concerning the effect of 
cone angle, B2O3 height, and ambient pressure indicate the importance of 
relatively high convective heat transfer at the crystal -ambient surface 
compared to the crystal -B2O3 surface. This finding disagrees with theoretical 
predictions of the relative heat transfer coefficients. The reconciliation of 
this discrepancy should be pursued to bring about a better understanding of 
the LEC crystal growth process and further reductions of the dislocation 
density. 
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Table I. Summary of EPO Measurements on LEG Ingots 


Ingot No. 

Wafer 

Location 

(IVinng— 

EPO icm~h 
(2) Center 

(3) Edge ' 

EPD Ratio 
(Center to Ring) 

1+ 

F 

.7.6 X 10^ 

4.6 X 10^ 

3.0 X 10® 

0.6' 


T 

6.1 X 1q5 

6.1 X 10® 

1,1 X 10® 

1.0 

2 

F 

2.2 X 10^ 

5.0 X 10^ 

2*3 X 10® 

2.3 


M 

3.2 X 10^ 

7.3 X 10^ 

2.5 X 10® 

2.3 

3 

F 

4.0 X 10^ 

8.0 X 10^ 

2.9 X 10® 

2.0 


T 

N/A 

N/A 

N/A 

- 


F 

4.0 X 10^ 

1.4 X lO^ 

4.0 X 10® 

3.5 


M 

3.0 X 10^ 

1.0 X 10^ 

N/A 

3.3 

5 

F 

1.5 X 10^ 

3.4 X 10^ 

1.7 X 10® 

2.3 


T 

1.2 X 1q5 

1.4 X 10® 

2.1 X 10® 

1.2 

6 

F 

1.8 X 10^ 

2.6 X 10^ 

8.0 X 10^ 

1.4 


T 

8.6 X 10^ 

7.7 X 10^ 

2.0 X 10® 

0.9 

7 

F 

l.O X 10^ 

2.5 X 10^ 

5.6 X 10^ 

2.5 


M 

2.5 X 10^ 

3.9 X 10^ 

7.8 X 10^ 

1.6 

8 

F 

1.4 X 10^ 

3.7 X 10^ 

1.0 X 10® 

2.6 


T 

N/A 

N/A 

N/A 

- 

9 

F 

1.4 X 10^ 

2.0 X 10^ 

2.5 X 10® 

1.4 


T 

1.0 X 105 

1.0 X 10® 

2.4 X 10® 

1.0 
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Table 

1 (Continued) 


Ingot No. 

Wafer 

Location 

^1) Ring 

EPD (cm'2) 
Jl) Center 

(3) Edge 

EPD Ratio 
(Center to Ring) 

10 

F 

1.0 X 10^ 

2.1 X 10^ 

1.1 X 10® 

1.9 . 


T 

N/A 

N/A 

N/A 

- 

11 

F 

7*5 X 103 

1.3 X 10^ 

1.9 X 10® 

1.7 


T 

8.1 X 10^ 

1.8 X lo5 

1.8 X 10® 

2.2 

12 

F 

1.2 X 10^ 

1.7 X 10^ 

2.5 X 10® 

1.4 


T 

9.0 X 10^ 

l.Q X 1 q5 

2.2 X 10® 

1.1 

13 

F 

3.5 X 10^ 

1.0 X 1q5 

1.5 X 10® 

2.9 


T* 

8.0 X 10^ 

1.1 X 10® 

2.0 X 10® 

1.1 

14 

F 

6.0 X 1q3 

1.8 X 10^ 

9.6 X 10^ 

3.0 


T 

N/A 

N/A 

N/A 

wm. 

15^^ 

F 

1.1 X 10^ 

2.4 X 10® 

2.7 X 10® 

2.2 


T 

1.3 X 105 

2.3 X 10® 

1.6 X 10® 

1.8 

16 

F 

1.3 X 10^ 

2.8 X 10^ 

1.7 X 10® 

2.2 


T 

1.4 X lO® 

2.2 X 1q5 

2.4 X 105 

1.6 

17 

F 

1.1 X 10^ 

2.0 X 10^ 

1.1 X lO® 

1.8 


T 

1.5 X 1q5 

2.5 X 10® 

1.7 X 105 

1.7 

18 

F 

8.5 X 1q3 

1.6 X 10^ 

1.2 X 10® 

1.9 


T 

9.7 X 10^ 

1.3 X 10® 

2.4 X 10® 

1.3 

♦: High EPD 

seed used 






+♦: FI at- top growth 
*: -3/4 ingot length area 
F: front, H: middle, T: tail 
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Table II. Effect of Cone Angle on Dislocation Density 


Ingot No. 

Cone Angle* 

— 

EPO (cm“2) 
Front 



1 

2.4 X 10® • 

15 

0 * 

2 

l.l X 10® 



3 

2.7 X 10® 



1 

4.0 X 10^ 

4 

10* 

2 

1.4 X 10^ 



3 

4.0 X 10^ 



1 

1.1 X 10^ 

17 

25* 

2 

2.0 X 10^ 



3 

1.1 X 10® 



1 

1.1 X 10^ 

10 

30* 

2 

2.1 X 10^ 



3 

in 

o 

K 

w«4 



1 

O 

w-4 

X 

GO 

■ 9- 

6 

50* 

2 

o 

X 

40 

• 

CM 



3 

8.0 X 10^ 



1 

1.4 X 10^ 

9 

62* 

2 

2.0 X 10^ 



3 

2.5 X 1q5 



1 

1.4 X 10^ 

8 

65* 

2 

3.7 X 10^ 



3 

1.0 X 10® 


*: Other growth parameters are similar. 


90 

J/3767A/cb 


> ♦ 


Rockwell Intematlonai 


OF POOR QUAUTY 


Table III. Effect of B 2 O 3 Height on Dislocation Density 



Height of*" 



EPD (cm*2) 


Ingot No. 

B 2 O 3 


Front 


Tan 



1 

3.5 X 

10 ^ 

8.0 X IG^* 

13 

270 gm 

2 

1.0 X 

10® 

1.1 X 10®* 



3 

1.5 X 

10® 

2.0 X 10®* 



1 

1.3 X 

10^ 

1.4 X 10® 

16 

390 gm 

2 

2.8 X 

10^ 

2.2 X 10® 



3 

1.7 X 

10® 

2.4 X 10® 



1 

1.2 X 

10^ 

9.0 X 10^ 

12 

500 gm 

2 

1.7 X 

10^ 

1.0 X 105 



3 

2.5 X 

10® 

2.2 X 10® 


: Other growth parameters are similar. 
: >^3/4 Ingot length area. 
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Table IV. Effect of Ambient Pressure on Dislocation Density 


Ingot No. 

Ambient* 

Pressure 


Front 

EPD (cm'2) 

Tall 



1 

6.0 X 

103 


14 

50 psi Ar 

2 

1.8 K 

10^ 

N/A'^ 



3 

9.6 X 

10^ 




1 

1.3 X 

10^ 

1.4 X lo' 

16 

300 psi Ar 

2 

2.8 X 

10^ 

2.2 X 10 



3 

1.7 X 

10® 

• 

X 

o 


*: Other sf'OMth parameters are similar. 
See text. 
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Effects of Seed Quality and Necking on Dislocation Density 


Ingot No. 

Necking* 

EPD (cm’2) 

Seed 

Tail 


No 


1 

7.6 X 10* 

1 

High 

2 

4.6 X 10^ 

• 


(5 X 10^) 

3 

3.0 X 10® 




1 

1.5 X 10^ 

5 

Yes 

High 

2 

o 

X 

o 

• 



(5 X IqS) 

3 

1.7 X 10® 

A 



1 

1.4 X 10^ 

9 

Yes 

Low 

2 

2.0 X 10^ 



(3.3 X lo3) 

3 

2.5 X 10® 

16 



1 

1.3 X 10^ 

No 

Low 

2 

2.8 X 10^ 



- 

(4.5 X 10^) 

3 

1.7 X 10® 

*: All cone angles >25* and other growth parameters are similar. 


C ' ^ 
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table VI. Effects of Melt Stoichiometry and Diameter 
Control on Dislocation Density 


Initial Melt* Diameter 
Ingot No. Composition Variation (mm) 


EPD (cm-2) 
Front Tall 


8 


10 


12 


16 


11 


S3.0t Ga 


51. St Ga 


51.5% Ga 


50.7% Ga 


Stoichiometric 


50.1% As 


50.3% As 


50.6% As 



1 

1.4 K 10^ 


± 4.0 

2 

3.7 K 10^ 

N/A 


3 

1.0 X 10® 



1 

1.8 X 10^ 

8.6 X 10^ 

t 3.0 

2 

2.6 X 10^ 

7.7 X 10^ 


3 

8.0 X 10^ 

2.2 X 105 


1 

1.4 X 10^ 

1.0 X 10® 

± 7.1 

2 

2.0 X 10^ 

1.0 X 105 


3 

2.5 X 10® 

2.4 X 10® 


0 

i 

1.1 X 10^ 


i 4.5 

2 

2.1 X 10^ 

N/A 


3 

1.1 X 10^ 



1 

1.5 X 10^ 

1.2 X 10^ 

t 8.5 

2 

o 

X 

o 

• 

1.4 X 10® 


3 

1.7 X lo5 

2.1 X 10® 


1 

1.2 X 10^ 

9.0 X 10^ 

t 1.6 

2 

1.7 X 10^ 

1.0 X 105 


3 

2.5 X 10® 

2.2 X 10^ 


1 

1.3 X 10^ 

1.4 X 105 

t 1.5 

2 

2.8 X 10^ 

2.2 X loS 


3 

1.7 X 10® 

2.4 X 105 


1 

7.5 X 103 

8.1 X 10^ 

t 1.5 

2 

1.3 X 10^ 

1.8 X 10® 


3 

1.9 X 10® 

1.8 X 10® 

growtn parameters are si mil ar. 
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FIGURE CAPTIONS 

Fig, 1 Slicing diagram for (100) LEC GaAs crystals. 

Fig. 2 Photograph of a KOH-etched (100) GaAs wafer showing the “ring" 

(Region 1), “center" (Region 2)c and "edge" (Region 3) regions, which 
have the lowest. Intermediate, and highest EPD, respectively. The 
dislocation density has fourfold symmetry Indicative of the <100> 
orientation. Also note the lineage polygonizatlon networks along 
<110>, which also have fourfold symmetry. 

Fig. 3 Photomicrograph of KOH-etched wafer quadrant showing microscopic View 
of dislocation density In regions 1-3. The EPDs In the edge and ring 
regions are lower along <110> than along <100>. 

Fig. 4 Radial dislocation density profiles across wafers obtained from the 

front, middle, and tall of a crystal. The radial profiles are 
"W"-shaped, and the average EPD Increases from the front to the tall. 

Fig. 5 Dislocation maps of longitudinal cross sections of cones with varying 

cone angles. 

Fig. 6 Dislocation maps of longitudinal cross sections of neck regions with 
varying neck-diameters: (a) 1.6 mn-dlameter neck showing severe 
deformation; (b) 3 mn-dlameter neck showing dramatic EPO reduction. 
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EFFECT OF MELT STOICHIOMETRY ON TWIN FORMATION IN LEC GaAs 


I 

I 

f! 

We have found that control over the melt stoichiometry is important I 

. , t: ■ 

to prevent twin formation in large-diameter, undoped, <10O> GaAs crystals 

i 

grown by the liquid encapsulated Czochralski (LEC) technique. Twenty GaAs 
crystals were grown from stoichiometric and non-stoichiometric melts. The 
incidence of twinning is reduced significantly when crystals are grown under 
As-rich, near-stoichiometric conditions. 

The undoped GaAs crystals used in this study were grown in a high 
pressure LEG growth apparatus. The pressure of the argon atmosphere during 

growth was 300 psi. The crystals were grown from pyrolytic boron nitride ' 

I 

(PBN) and quartz crucibles. The growth direction was <100>. The water | 

content of the B 2 O 3 encapsulant (Johnson Matthey) material ranged from 160 to 

500 ppm. The weight and diameter of the crystals were typically about 2.4 Kg 

and 3 in., respectively. The crystals were grown with cone angles varying 

from 0 to 70*. The cone angle is defined as the angle between the wall of the 

cone and the horizontal; a crystal with a 0 * cone angle is flat-topped. 

I ' ■ 

'h ' , 

' I 

Kay Worde: liquid enoapaulated Cecsahralaki growth, twinning, malt atoiohio^ | 

matry, gallim areanida, | 

I ■ 
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The GaAs meTt was prepared by in-situ synthesis ( 1 ) starting with a 
3 Kg charge of high purity As (Cominco) and Ga (Kawecki). The stoichiometry 
of each melt was changed by varying the composition of the initial charge. To 
make an accurate determination of the initial melt composition for crystals 
grown from PBN crucibles, it was necessary to take into account the loss of As 
from the charge during the heat-up cycle. We found that the loss of arsenic, 
which varied from about 10 to 90 g, resulted from incomplete wetting of the 
B 2 O 3 to the PBN crucible before synthesis. The weight loss was determined by 
comparing the weight of the initial charge with the weight of the crystal and 
the charge remaining in the crucible after growth. It is important to note 
that the B 2 O 3 was separated from the remaining charge before determining the 
weight loss. This technique eliminated possible errors associated with the- 
evaporation of B 2 O 3 during growth. We emphasize that the presence of excess 
Ga(or As) in the remaining charge always correlated with the melt stoichiometry 
calculated on the basis of the weight loss measurements. The composition of 
the melt at the beginning of growth was effectively varied from 0.46 to 0.51 As 
atom fraction. The melt composition at the termination of growth was 
determined by adjusting the initial melt composition for the crystal weight. 

We did not determine the stoichiometry of melts synthesized in quartz 
crucibles. However, in subsequent experiments we found that only from 5 to 20 
g of As is typically lost from a stoichiometric 3 Kg charge contained in a 
quartz crucible during heat-up. Again, B 2 O 3 was separated from the GaAs 
residual charge for these weight loss measurements. In addition, when quartz 
crucibles were used, we observed relatively small amounts of As deposited on 
the walls of the growth chamber compared to growth with PBN crucibles, and no 
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excess Ga in the remaining charge. Since stoichiometric charges were always 
used with quartz crucibles in the present study, we took the initial melt 
stoichiometry to be slightly Ga>rich with an estimated composition of 0.498 As 
atom fraction. 

The results of our study, summarized in Table 1, show that the incidence 
of twinning is significantly reduced when crystals are grown from As-rich melts. 
Only 4 of 12 (33%) crystals grown from Ga-rich melts were single. On the other 
hand, 7 of 8 ( 88 %) crystals grown from As*rich melts were single. Furthermore, 
the incidence of twinning could not be correlated with other growth parameters, 
such as the wetness of the 8203 ( 1 ), the cone angle (see Table 1), or fluctu- 
ations in the diameter of the crystal. The results indicate a sharp increase 
in t«^if?ning probability on the Ga-rich side of the stoichiometric composition. 

Previous studies (2,3) have shown that the incidence of twinning in 
small -diameter GaAs crystals can be reduced by growing with gradual cones; i.e., 
large cone angles. Our results show no correlation between the incidence of 
twinning and cone angle in large -diameter crystals. Moreover, the significantly 
reduced incidence of twin formation experienced using As-rich melts in the 
present study was achieved with small com angles ranging from 0 to 35®. 

Growth experiments employing quartz crucibles were not conducted with 
As-rich, undoped melts to compare with the results obtained with the Ga-rich 
melts. However, we have recently grown several crystals doped with Se, Si, and 
Zn from As-rich melts using quartz crucibles. The incidence of twin formation 
was very low in this series of experiments (8 out of 9 crystals were sing’n). 

Our combined results therefore indicate that twin formation is independent of 
the crucible material. 
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The twinned crystals were categorized according to the twin 
morphology. One group was characterized as having only one longitudinal twin, 
which nucleated at the surface of the crystal and cut the crystal obliquely on a 
(111) plane. The twinned region of one such crystal was found by x-ray analysis 
(4) to be oriented with the <112> direction parallel to the growth direction. 

The second group of crystals contained multiple twins. Twins in all crystals 
Invariably nucleated at one of the four peripheral facets that run axially along the 
crystals. (The peripheral facets result from the Intersection of (111) As and 
(111) Ga facet planes with the edge of the crystal along <110> directions that 
are perpendicular to the <10G> growth axis). No preference was observed for 
either As or Ga peripheral facets as nucleation sites for twins* 

The reduced Incidence of twin formation in As-rich melts has been 
reported for GaAs grown by the Bridgman (5) and modified Gremmelmaier (2) tech- 
niques. The consistent effect of melt stoichiometry on twin formation in GaAs 
grown by three different techniques may therefore reflect fundamental behavior 
of the material. The dramatic variation in the incidence of twinning over a 
relatively small range of melt compositions observed in the present study 
suggests that the stoichiometry of the solid at the growth interface could play 
an important role. Thus, the variable resistance of the crystal to twin 
formation could be related to different solidification kinetics depending on 
whether vacancies. Interstitials, or antisite defects incorporate Into the 
solid. 

In summary, we have shown that the incidence of twin .formation in 
large-diameter, undoped, <100> LEG GaAs Is reduced when the melt composition is 
slightly As-rich. In view of the potential for the loss of As from the charge 
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when using in-situ synthesis, the yield of single <100> crystals will depend on 
close control of the melt composition. Final ly. the results suggest that the 
barrier to twin formation is related to the stoichiometry of the solid at the 
solidification front, 
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